Entropic Fluctuations in Statistical Mechanics I. 
Classical Dynamical Systems 



o 

o 

Q ! V. Jaksic 1 , C.-A. Pillet 2 , L. Rey-Bellet 3 

O 

t— i ■ 

1 Department of Mathematics and Statistics 
McGill University 
805 Sherbrooke Street West 
■ Montreal, QC, H3A 2K6, Canada 

2 Centre de Physique Theorique* 

Universite du Sud Toulon- Var, B.P. 20132 

CN ■ F-83957 La Garde Cedex, France 

> 

■"sj - , 3 Department of Mathematics and Statistics 

Lederle Graduate Research Tower, Box 34515 
University of Massachusetts 
Amherst, MA 01003-4515, USA 

O ' 



December 13, 2010 
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statistical mechanics. Our main objective is to display the minimal, model independent mathematical structure 
at work behind fluctuation theorems. Besides its conceptual simplicity, another advantage of our approach is 
its natural extension to quantum statistical mechanics which will be presented in a companion paper. We shall 
discuss several examples including thermostated systems, open Hamiltonian systems, chaotic homeomorphisms of 
compact metric spaces and Anosov diffeomorphisms. 
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1 Introduction 

This is the first in a series of papers devoted to the so-called Fluctuation Theorems of non-equilibrium statistical 
mechanics. This series is a part of the research program initiated in [Pi JP1 JP2| that concerns the development 
of a mathematical theory of non-equilibrium statistical mechanics within the framework of dynamical systems. 

The first fluctuation theorem in statistical mechanics goes back to 1905 and the celebrated work of Einstein on 
Brownian motion. The subsequent historical developments are reviewed in |RM| (see also the monographs |GM 
IKTH1 ) and we mention here only the classical results of Onsager MOnll lOn2 l. Green IGrlllGr2l . and Kubo IKubl 
which will be re-visited in this paper. Virtually all classical works on the subject concern the so called close to 
equilibrium regime in which the mechanical and thermodynamical forces (affinities) are weak. One of the key 
features of modern fluctuation theorems, suggested by numerical experiments |ECM| and established theoretically 
for the first time by Evans and Searles OESI and by Gallavotti and Cohen 0GClllGC2l . is that they hold for systems 
arbitrarily far from equilibrium and reduce to Green-Kubo formulas and Onsager relations in the linear regime near 
equilibrium. The seminal papers [ECM ES GC1 GC2 1 were followed by a vast body of theoretical, numerical and 
experimental works which are reviewed in |RM|. The Evans-Searles (ES) and Gallavotti-Cohen (GC) Fluctuation 
Theorems are the main topics of the present work. 

The basic two paradigms for deterministic (dynamical system) non-equilibrium statistical mechanics are the so 
called thermostated systems and open systems. Thermostated systems are Hamiltonian systems (with finitely many 
degrees of freedom) driven out of equilibrium by an external (non-Hamiltonian) force and constrained by a de- 
terministic thermostating force to stay on a surface of constant energy. Open systems are Hamiltonian systems 
consisting of a "small" Hamiltonian system (with finitely many degrees of freedom) interacting with, say two, 
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"large" reservoirs which are infinitely extended Hamiltonian systems. The reservoirs are initially in thermal equi- 
librium at distinct temperatures and the temperature differential leads to a steady heat flux from the hotter to the 
colder reservoir across the small system. Throughout the main body of the paper we shall illustrate our results on 
an example of thermostated system and an example of open system. 

The majority of works on fluctuation theorems concern classical physics. In the quantum case comparatively little 
is known and there are very few mathematically rigorous works on the subject (see IITMI IDDM1 IRol |Ku21). The 
present paper, which concerns only the classical case, originates in our attempts to find a proper mathematical 
framework for the extensions of ES and GC Fluctuation Theorems to quantum physics. One of the difficulties in 
finding such a framework stems from the fact that it was already lacking at the classical level. Indeed, even the 
basic examples of thermostated systems and open systems were studied in the literature in an unrelated way and it 
was far from obvious which aspects of the theory are model dependent and which are universal. For example there 
was no clear universal rationale in the choice of the "entropy production" observable (also called "phase space 
contraction rate" or "dissipation function") which plays a central role in the theory. 

A model independent definition of the entropy production has been proposed by Maes, in the context of stochastic 
(Markovian or Gibbsian) dynamical systems, see e.g. BMa2l . We take here a different and complementary route 
and discuss non-equilibrium statistical mechanics within the context of deterministic dynamical systems. Our work 
is mostly of a review nature and we do not prove any new specific results. Rather we organize the existing set of 
ideas and results in an axiomatic abstract framework that unifies virtually all deterministic models discussed in 
the literature (in particular, open infinite systems and thermostated finite dimensional systems will be treated in a 
unified manner) and clarifies the mathematical structure of the theory. The framework has a direct extension to 
non-commutative dynamical systems and in particular to quantum mechanics and this will be the subject of the 
remaining papers in the series. Our principal new results concern the quantum case and we will focus here only on 
those aspects of the classical theory that can be extended, within the framework of dynamical systems, to quantum 
statistical mechanics. 

We have made an attempt to expose the results in a pedagogical way and the only prerequisite for the principal part 
of the paper is a basic knowledge of probability and measure theory. 

The paper is organized as follows. 

In Section |2] we introduce our dynamical system setup and review the properties of relative entropies that we will 
need. In Section [3] we introduce the basic objects of the theory, the entropy cocycle and the entropy production 
observable, discuss their properties, and prove the finite time Evans-Searles Fluctuation Theorem. The results 
described in this section hold under minimal regularity assumptions that are satisfied in virtually all models of 
interest. 

In Section [4] we start the discussion of thermodynamics by introducing control parameters (mechanical or ther- 
modynamical forces) to our dynamical system setup. The finite time Evans-Searles Fluctuation Theorem is then 
generalized to this setting. Following the ideas of Gallavotti MGall (see also |LS2 1) we use this generalization to 
derive finite time Green-Kubo formulas and Onsager reciprocity relations. 

The results of Sections [3] and |4] concern the system evolved over a finite interval of time and are very general. In 
particular, they do not require any ergodicity assumptions. Section|5]concerns the large time limit t — » oo. Under 
suitable ergodicity assumptions we derive Evans-Searles Fluctuation Theorem on the basis of its universally valid 
finite time counterpart and prove the Green-Kubo formula and Onsager reciprocity relations. 

Section [6] is devoted to the Gallavotti-Cohen Fluctuation Theorem. After introducing the key concept of non- 
equilibrium steady states (NESS), the GC Fluctuation Theorem is stated as, essentially, an ergodic-type hypothesis 
concerning the NESS and the entropy production observable. The Green-Kubo formula and Onsager reciprocity 
relations also follow from the GC Fluctuation Theorem. 

One advantage of our abstract axiomatic framework is that it allows for a transparent comparison between the 
ES and GC Fluctuation Theorems. It turns out that from the mathematical point of view these two theorems are 
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equivalent up to an exchange of limits (see Relation j7.62\ ). This exchange of limits may fail even in some very 
simple models and its validity can be interpreted as an ergodic property of the underlying dynamical system. We 
raise this point to the Principle of Regular Entropic Fluctuations which is introduced and discussed in Section|7] 
After this work was completed we have learned that related ideas have been previously discussed in [RM|. 

Sections IMTTI are devoted to examples. In Section|8]we discuss several toy models which illustrate the optimality 
of our assumptions. In Section [9] we develop the non-equilibrium statistical mechanics of Gaussian dynamical 
systems. Chaotic homeomorphisms of compact metric spaces are discussed in Section [10] Finally, in Section [TT] 
we discuss the non-equilibrium statistical mechanics of Anosov diffeomorphisms of compact manifolds. In each 
of these examples we verify the validity of the proposed principle of regular entropic fluctuations. 

For the convenience of the reader, a table of frequently used abbreviations and symbols is provided on page [67] 

The ergodic-type hypotheses introduced in this paper typically concern existence of certain limits as time t — > oo, 
the regularity (differentiability, etc) properties of limiting functions w.r.t. control parameters, and the validity of 
exchange of order of limits and derivatives. The introduced hypotheses are minimal (i.e., sufficient and necessary) 
to derive fluctuation theorems and their implications from the universally valid structural theory discussed in Sec- 
tions 2-4. The verification of these hypotheses in concrete models leads to a novel class of (analytically difficult) 
problems in ergodic theory of dynamical systems. 

Acknowledgment. The research of V.J. was partly supported by NSERC. The research of C.-A.P. was partly 
supported by ANR (grant 09-BLAN-0098). The research of L.R.-B. was partly supported by NSF We wish to 
thank C. Liverani and D. Ruelle for useful discussions. The results of this paper were presented by its authors in 
mini-courses at University of Cergy-Pontoise, Erwin Schrtidinger Institute (Vienna), Centre de Physique Theorique 
(Marseille and Toulon), University of British Columbia (Vancouver), Ecole Polytechnique (Paris), Institut Henri 
Poincare (Paris) and Ecole de Physique des Houches. The paper has gained a lot from these presentations and we 
wish to thank the respective institutions and F. Germinet, J. Yngvanson, R. Froese, S. Kuksin, G. Stoltz, J. Frohlich 
for making these mini-courses possible. 



2 Basic notions 

2.1 Phase space, observables, states 

Let M be a set and T a cr-algebra in M. We shall refer to the measure space (M, T) as the phase space. If M is a 
topological space, we shall always take for T the Borel cr-algebra in M. 

An observable is a measurable function / : M —¥ C and we denote by 0(M) the complex vector space of all ob- 
servables. B(M) denotes the subspace of all bounded observables. Together with the norm ||/|| = sup xeM \f(x)\, 
B(M) is a Banach space. If M is a topological space, C(M) denotes the Banach space of all bounded continuous 
observables. The corresponding spaces of real valued observables are denoted by Om(M), Br(M), Cr(M). 

A state is a probability measure on (M, J 7 ) and S denotes the set of all states. The expectation value of an 
observable / w.r.t. the state v is denoted by 



If f = (/i, . . . , /at) is a vector-valued observable we set u(f) = (f(/i), ■ ■ ■ , ^(/jv)). We shall equip S with the 
weakest topology w.r.t. which the functionals 59f4 v(f) are continuous for all / G B(M). If 6 : M —> M is 
a measurable map, we denote by v o O^ 1 the measure T 3 A M> v(9^ 1 (A)). Clearly v o 6~ 1 (f) = v(f o 6) for 
all / e B(M). A map 9 : M —> M is called involutive if d o 9(x) = x for all x G M. 
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We shall say that a state v is normal w.r.t. uj G S iff v is absolutely continuous w.r.t. uj (denoted v <C oj). The set of 
all states which are normal w.r.t. ui is denoted j\f u . Two states v and u are called equivalent iff v <C w and cj <C ^, 
i.e., iff v and w have the same sets of measure zero. The Radon-Nikodym derivative dv/du), which will play an 
important role in this paper, is defined as an element of L 1 (M , duo) and is an equivalence class of functions rather 
than a single function. For this reason the notion of observable is extended as follows. Given a state uj, let = 
{/ € O(M) | f(x) = for oj-a.e. x} and let 0(M) U = 0(M)/Z u be the quotient vector space (the elements of 
0(M) U are equivalence classes w.r.t. the relation f^g-^f — gE Z u ). Similarly, L°°(M, duj) = B(M)/Z U . 
As usual in measure theory, dealing with equivalence classes instead of single functions is natural and causes no 
difficulties, the classes are called functions, etc. 



In what follows, we adopt the shorthands 



dv 



2.2 Relative entropies 

The relative entropy of a state v w.r.t. a state to is defined by 

n i/ I \ J -oo if u#M,, 

Ent(y\u) = < (2.1) 
[-"(tu\u) rfueAfu,- 

Since — t v \ u < A~^ — 1 and u(A~}) — 1, relative entropy is well defined as a map from S x S to [— oo, 0]. Its 
basic properties are (see e.g. [OP]): 



Theorem 2.1 (1) Foruj, v e 5, Ent(^|w) = inf/ 6 s H (M) [l°gw(e-' ? ) — v{f)\. 

(2) Forwe5 and / G B R (M), logw(e^) = sup„ e5 [Ent(V|w) + z/(/)]. 

(3) Concavity: For wi, ti>2> ^li ^2 G 5 ana! A G [0,1], 

Ent(A^i + (1 - A)iA,|Acji + (1 - \)uj 2 ) > AEnt(i/i|wi) + (1 - A)Ent(i/ 2 |w 2 ). 

(4) Ent(i/|w) < 0/or all uj, v G 5 ant/ Ent(j/|w) = if and only if v = uj. 

(5) : M —¥ M is a measurable bijection, then ~Ent(i/ o ^ 1 |oj o = Ent(^|w). 

(6) The relative entropy is an upper semicontinuous map from S x S to [—oo, 0], that is 

Ent(i/|w) > limsupEnt(i/ Q |w Q ), 

a 

/or aW convergent nets v a v and UJ a — > uj in S. 

(7) For any uj G S and any finite constant C, the set \v G S | Ent(z/|w) > C} is compact in S. 

The Renyi relative entropy of order a G R, BReL is defined by 

Ent a (i/w)=^ /A N * r (2.2) 
" V 1 ' \log W (A«J if^GA^. 



This generalization of relative entropy has found numerous applications (see |BS OP] for references and additional 
information). We list below several properties of the Renyi relative entropy that are relevant for our purposes: 
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Proposition 2.2 Suppose that v £ J\f u . 

(1) If 9 : M — > M is a measurable bisection, then Ent Q (i^ o Q~ x \io o = Ent a (p\u)). 

(2) R 3 a i— > Ent a (i/|w) G] — oo, oo] is a convex function. It is real analytic and non-positive on ]0, 1[. It is 
positive for a [0, 1]. 

(3) lim — — Ent a (i>\u) = Ent(v\uj). 
atl 1 — a 

In the remaining statements we assume that v and uj are equivalent. 

(4) EntoH^) = Enti(i/|u>) = 0. 

(5) ~EiTti a (i/\ui) = Enti_ Q (u;|i/). 

(6) Ent a (i/|w) > aEnt(w|^). 

2.3 Dynamics 

Let I be an index set whose elements are interpreted as instances of time. We shall always assume that X = Z 
(the discrete time case) or I — R (the continuous time case). A dynamics (j> = {0* 1 1 £ 1} on M is a group of 
invertible measurable transformations <j? : M M describing the evolution of the system. More precisely, we 
shall assume: 

(Fl) 4>° is the identity map and <p t+s = 0* o <f> s for all s, t € I. In particular, for all t E I, (ft is an 
automorphism of the measurable space (M, F). 

(F2) The map (t, x) h- > (^(x) is measurable. 

The assumption (F2) is relevant only in the case I = R (in this case, the dynamics <p is a flow on M ). In the 
discrete time case I = Z, the dynamics is obtained by iterating the time 1 map <f> = (f) 1 and its inverse <f> . We 
will sometimes write <fi n instead of 0*. 

A dynamics <j) on M induces transformation groups on O(M) andSby f t — /o0*, v t — vo$r l . They are clearly 
related by vt{f) — v {ft)- A state v is called steady (or stationary) if v t — v for all t. We denote by Si the set of 
all steady states. 

2.4 Reference state 

The starting point of our discussion is a classical dynamical system (M , <fi, uj), where is a given dynamics on M 
and oj a given reference state satisfying the following regularity assumption: 

(C) ujt and cj are equivalent for all tel. 

In non-trivial models that arise in non-equilibrium statistical mechanics J\f u OSi =0. In particular, lu Si. In this 
important aspect our starting point differs from the usual one in the ergodic theory of classical dynamical systems 
where the reference state cj is assumed to be invariant under the dynamics. 

Assumption |(C)| ensures that <f> preserves and hence naturally induces a group of transformations of C(Af) w 
andi°°(M,dw). 

Assumptions |(F1)| |(F2)| and |(C)| are our fundamental working hypothesis and will be assumed in the following 
without further notice. 
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2.5 Time reversal invariance 

A time reversal of the dynamics <f> on M is an involutive measurable transformation i? : M — > M such that 

i? o ft = (/>-* o i? for all t G I. 

A state w G 5 is called time reversal invariant (TRI) if u o i? = oj. In this case preserves Z w and induces an 
involution on 0{M) U and L°°(M, da;). Note that if w is TRI, then uj t o ?9 = w_ t . 

The dynamical system [M, <j>, u) is called TRI if M is equipped with a time reversal i? of such that w is TRI. 

Time-reversal invariance will play a central role in our discussion. Other symmetries can have important conse- 
quences on statistical properties of the dynamics, e.g. the conformally symplectic structure of some systems leads 
to symmetries in their Lyapunov spectrum (see MWL1 IMDI). Such symmetries, however, will not play a role in our 
work. 

If the system (M, <fi, w) is not TRI, for the purpose of model building the following construction is useful. Set 

M = MxM, $(x,y) = ((j) t (x),4>' t (y)), dw = dto <g> dw. 
Then (M, <f>, ui) is TRI with the time reversal $(x, y) = (y, x). 

3 Finite time entropy production 

3.1 Entropy cocycle 

Since ui t and cj are equivalent measures, 

is well-defined. It satisfies the following additive cocycle property: 
Proposition 3.1 For all t, s G I one has 

c t+s = c s + c* o 4> s . 

In particular c° — and c~ l = — c* o </> _t . 

Proof. We adopt the shorthand A* = A Ut | w . For / G B(M) and s,t el one has uj t+s (f) = uj(A t+s f) and 
«*+.(/) - ".(/*) = ^(A s ft) = u((A s o 0-' /) o 0') = Wt (A s o 4r* /) = W (A S o A* /). 

Hence, 

A t+S = A s of A*, (3.3) 
where the equality is in 0(M) U . Taking the logarithm we derive 

Our first identity follows immediately. The second one follows from the substitution s = and the third one is 
obtained by setting s = —t. □ 

We shall call c* the entropy cocycle of the dynamical system (M, </>, w). The entropy cocycle of a lTRII dvnamical 
system enjoys the following additional property. 
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Proposition 3.2 If the system (M, <p, uj) is YYKU with a time reversal i9, then 

c t od = c - t , (3.4) 

holds for all iel. 

Proof. Setting again A* = A Ut | u we have, for any / G B(M) and tel, 

w(A'otf /) = wotf (A* /otf) = w(A*/°i?) = w t (/oi?) = w(/o#o</>*) = w (/o^-*oi9) = w(/o0-*) = u>(A -t /)- 

The resulting identity A - * = A* o $ further leads to A~* o <^> _ * = A* o $ o _t = A* o <^>* o Taking the logarithm 
gives the result. □ 

3.2 Entropy balance equation 

By Definition ( 12.1b of the relative entropy one has 

Ent(wt|w) = -wt(4; t | w ) = -w(4, t | w ° 0') = -w(c*). 
Since Ent(o;|a;) = this identity can be rewritten as 

w (E*) = -- (Ent(cj 4 |w) - Ent(w|w)) , 

where 




(3.5) 



Thus, we can interpret E* as the observable of mean entropy production rate over the time interval [0, t]. We shall 
call the relation 

Ent(w t |w) = -tw(E*), (3.6) 
the entropy balance equation. Its immediate consequence is the important inequality 

w(E*) G [0,00], (3.7) 

which holds for all i > 0. 
The cocycle property yields 




We note for later reference that if (M, 4>, ui) is ITR1I then Proposition 13 . 21 further leads to 

E* o d = — = -E"* = -E* o cjr*. (3.8) 

3.3 Finite time Evans-Searles symmetry 

Let P t be the law of the real-valued random variable E t , i.e., the Borel probability measure on R such that for any 

/ e B(R), 

P t (/)= W (/(S*)). 

Let t : R — > R be the reflection t(s) = — s and define the reflected measure P t = P t o r. 
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Proposition 3.3 If (M, <j), uj) is \TRl\ then, for any t G X, the measures P t and Pt are equivalent and 

gw=e- ,3.9, 



Proof. For / e B(R), Equ. (O and the fact that uj t o ■& = uj_ t yield 

Pt(f) = "(/(-£')) = <*(/(-£* o 0-')) = o 0)) = «_ t (/(E*)) - c(e-* s V(E 4 )), 

and the statement follows. □ 

To our knowledge, the relation ( 13. 9t was first obtained by Evans and Searles in [ES| and is sometimes called 
the transient fluctuation theorem. We shall call it the finite time ES-identity. We stress its universal character: 
besides the ITRII assumption it only relies on the minimal hypothesis |(F1)| |(F2)| and |(C)| In a loose sense, it 
can be understood as a dynamical form of the second law of thermodynamics: on the finite time interval [0,t], 
the probability to observe a negative mean entropy production rate — s is exponentially small compared to the 
probability to observe the positive value s. 

The ES-identity can be re-formulated in terms of Renyi entropy. For a e I, we adopt the shorthand 

e t (a) = Ent Q (a; t |w) = logw(e Q£ "ti-) = loga;(e Qts_t ). (3.10) 
By TheoremO(2), if £ Ut \ u G L°°(M, dcj), then 

e t (a) = sup [Ent(z/|w) + av(l Ut \ u )]. (3.11) 

This variational characterization will play an important role in the extension of the theory of entropic fluctuations 
to non-commutative dynamical systems. 

The basic properties of the functional ( 13.101 ) follow directly from Proposition ^. 21 We list them for later reference: 
Proposition 3.4 (1) For alltEl the function 

l9a4 e*(a) G] — oo, oo], 
is convex, satisfies et(0) = et(l) = and 

j e t (a) e] - oo,0] j/aG[0, 1], 
I e t (a) G [0, oo] otherwise. 

(2) It satisfies the lower bound 

e t {a) > min(aEnt(a;|a; t ), (1 - a)Ent(w t |a;)) . (3.12) 

(3) It is real analytic on the interval ]0, 1[. 

Remark. The relation e t (l) = is sometimes called the non-equilibrium partition identity or Kawasaki identity, 
see IICWWI . 

Note that if the system is lTRTl Equ. (|3T8) implies w(e Qts_t ) = w(e" Q * sfo,? ) = w(e" atst ) so that 

e t (a)=log W (e- Qtst ) . 
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Proposition 3.5 (1) For any t S land a e M. one has e t (a) — e_t(l — a). 
(2) If(M,(j>,ui) is \TRl\ then e_ t (a) = e t (a) and hence 

e t (a) = e t (l - a). (3.13) 

Proof. Parts (1) and (5) of Proposition ^. 2l implv 

et(a) = Ent Q (w t |w) = Enti_ Q (w|w t ) = En.ti_ a (a;_f|a;) = e_ t (l - a). 
Since lTRll implies w t oi? = a>_ t , Part (1) of Proposition ^. 2l allows us to conclude 

e_t(a) = Ent Q (w_t|w) = Ent a (cJt o $|cj o ■&) = Ent a (a;t|a;) = et(a). 

□ 

We shall call Relation ( 13.13b the finite time ES- symmetry. We finish this section with the observation that the finite 
time ES-symmetry is an equivalent formulation of the finite time ES-identity. 

Proposition 3.6 For each t G X, the following statements are equivalent: 



(1) The measures Pt and Pt are equivalent and satisfy the ES-identity \3. 9i . 

(2) For all a G R, e_ t (a) = e_ t (l - a). 

Proof. It suffices to notice the relation between the functional e_t(a) and the Laplace transform of the measure 
Pt . One has 

e_ t (a) = log^ (e" ats *) - log J e ats dP t (s), 

and hence 

e_*(l - a) - log J e-^ ts dP t ( S ) = log J e ats e- ts dP t (s). 

□ 



3.4 Entropy production observable 



For a discrete time dynamical system the cocycle property c + =c|c of implies 



s=0 

where 

<x = c 1 =l wAu oci>. (3.15) 
In particular, we can express the mean entropy production rate observable as 



t t ^ 

s=0 
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Consequently, the entropy balance equation d3.6t becomes 



t-i 



Ent(wt|w) = — u)(a s ). 



s=0 

We shall call a the entropy production observable. For obvious reasons the entropy production observable is also 
often called phase space contraction rate. 

Basic properties of the entropy production observable are: 

Proposition 3.7 uj(ct) > 0, w_i(cr) < 0, and if (M. 6. u) is \L'Rl\ then got} = — g_i . 
Proof. —u>(a) = —<^i(i Ul \ui} = Ent(cL>i|w) < implies to(a) > 0. Jensen's inequality 

e"- 1 M <U- 1 (e")= ( j(A Ul{u ) = l i 
implies oj_i(ct) < 0. The last statement follows from d3.4t and the cocycle property r 1 = -c 1 o D 

It is not possible to define the entropy production observable of a continuous time dynamical system at the current 
level of generality. We shall make some minimal regularity assumptions to ensure that the entropy cocycle has a 
generator a, i.e., that the continuous time analog of Equ. d3.14| i holds. 

(El) 1. The function R 3 1 1-> A Wt[w € L l {M, du) is strongly C 1 . 
2. The entropy production observable 



= a; A 



ui t \ui 



t=0 

is such that the function R 3 t H> cr t € L 1 {M 1 du) is strongly continuous. 

Remark. If M is a complete, separable and metrizable space, then the Koopman operators T* : / i-» A U( | u /_( 
form a strongly continuous group of isometries of L 1 (Af, dw). Denote by i its generator, i.e., T* = e tL . Since 
A Wt | u — T l l, Part 1 of Assumption |(El)| is equivalent to 1 G Dom(L), the domain of L, and then <r = LI. 

Proposition 3.8 Suppose that (El) holds. Then: 

(1) For all x G M the function 1 1— > A LJt | w (a;) is absolutely continuous and 

holds for uj-almost all x G M ant/ Lebesgue almost all t G R. 

(2) For aZZ i G R f/ze identities 

A i = eJo CT - sCU 

anc/ 

„t 



c ( = / a s ds, (3.16) 

yrfin 0(M) W . 
(3) w(ct) = 0. 
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(4) If(M, cf>, u) is HM then ao-d = -a. 
Proof. We again set A* = A Wt i w . 

(1) The cocycle property d3.3l ) and Part (1) of Assumption |(El)| yield 

w(|A t+s - A* - sAV_ t |) = ui (A* |A S o (j)- 1 - 1 - so-_ t |) = cj(|A s - 1 - sa\) = o(s), 
as s — > 0, from which we conclude that 

d_ 
~dl~ 

holds strongly in L 1 (M, dw). By Part (1) of Assumption |(E 1 )| we can assume that for any x <E M the function 
1 1— > A* (x) is absolutely continuous and that for w-almost all x E M 

^A t (x)=A t (x)a. t (x), 

holds for Lebesgue almost all i e R (see e.g. Theorem 3.4.2 in [HP|). 

(2) By Part (2) of Assumption |(El)| the Riemann integral 

It = / cr-s ds, 
Jo 

defines a strongly C 1 function t >-> G L 1 (M, dw). As before, we can assume that £ i-> £j(a;) is absolutely 
continuous for all x £ M and that for w-almost all x € M 

^t(x) = cr- t (x), 

holds for Lebesgue almost all i 6 1. Consequently, for all x € M the function t H> F t (x) = A t (x)e~ it ^ is 
absolutely continuous and for w-almost all x € M 

^1(*)=0, 

holds for Lebesgue almost all t G R. We conclude that for alU € R, i*t(a:) = F (x) = 1, i.e., A*(x) = e^ t(;r ) for 
w-almost all x G M. Equ. ( 13.16b follows immediately. 

(3) Differentiating the identity e f (1) = w.r.t. t at i = we derive w(<r) = 0. 

(4) Follows from the identity d3~4l . □ 

Remark. Under the strong continuity condition of Assumption |(E1)| the identity ( 13.16b holds in 0(M) U with a 
Lebesgue integral. It also holds in L 1 (M, dw) with a strong Riemann integral. 



The relation ( 13.16b yields 



1 f* 
S+ = — / (T s ds, 



and so the entropy balance equation and the ES-functional can be written as 



Ent(tJf |cj) = — / w(<7 s )ds, 
Jo 
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e t (a) = logw (e% t<T - 3 
For lTRll svstems one has <j_ s o i? = — <t s and in this case 

e t (a) = logo; (V tt J>= d ^ . 

Unless otherwise stated, in the sequel we will only consider continuous time dynamical systems. The discrete time 
case is very similar, time-integrals being replaced by appropriate sums. 

3.5 L p -Liouvilleans 

In this section, in order to avoid unessential technicalities, we shall assume in addition to |(El)[ 
(E2)cr € L°° (M, duj). 

Let p e] - oo, oo], p 7^ 0, and / £ L°°(M, doj). We shall consider the following special class of Ruelle transfer 
operators 

U p (t)f = A* ,„/_, = e7ti°-° ds f. t . (3.17) 

One easily shows: 

Proposition 3.9 Under Assumptions |(E1)|(E2)] Equ. ( 15.771 ) defines a family of bounded linear operators on 
L°°(M, duj) which satisfies: 

(1) U p {0) = I andU p (t + s) = U p {t)U p (s). 

(2) Ifp- 1 + q- 1 = 1, then u{[U p {t)f][U q {t)g]) = u(fg), 

(3) Ld(\U p (t)f\ p ) — uj(\f\ p ). For p £ [l,oo], U p (t) extends to a group of isometries of L p (M,du)). 

(4) U p (t) extends to a group of bounded operators on L 2 (M, du) such that U*{t) = U q (—t) and 

||C/ p (t)||<el'l^, m p = ||a|U. 

(5) Suppose that U p (t) is strongly continuous on L 2 (M,dui) and let L p be its generator, U p {t) = e tLp . Then 
L* = —L q , sp(Lp) C {z | |Rez| < m p }, Dom(L p ) = Dom(i oc ), and for f £ Dom (L p ) 

Lpf = Loo/ + -/■ 
P 

We shall call the operator L p the L p -Liouvillean. In Part (5), Dom(A) denotes the domain of the operator A and 
sp(^4) its spectrum. 

If a = 1/p, then 

e t (a) = log(l,e tL n) = log / e'Mdu;. (3.18) 

J M 

Similarly to ( 13.111 ). this operator characterization of e t (a) will play an important role in the extension of the theory 
of entropic fluctuations to the non-commutative setting. 



15 



If a is unbounded, the operator U p (t) could be unbounded and Proposition [3]9] may not hold as formulated. The 
technical aspects of its extension are then best carried out in the context of concrete models. 

In the discrete time case the concept of L p -Liouvillean is not very natural and instead one deals directly with the 
transfer operator 

U p f = e^f_ 1 . (3.19) 
If a G L°°(M, dw), then Parts (l)-(4) of Proposition |3~9l obviously hold, U p (n) = U% and e„(a) = log(l, U%1). 

3.6 Examples: Differentiable dynamics and thermostated systems 

Let U C R n be an open connected set and let cf> : U — » U be a C 1 -diffeomorphism. Denote by D<f> its derivative. 
Let M C U be compact and suppose that 4>{M) C M. Finally denote by cj the normalized Lebesgue measure on 
M. The entropy production observable of the discrete time dynamical system (M, tj>, w) is given by 



a = -log|detD(/)| 



To describe the continuous time case let X be a C 1 -vector field on U. Assume that the flow <f> generated by the 
differential equation 

—x* = X{x t ). 
dt y th 

satisfies 4>\M) C M. Then <j) is C 1 on R x AI. The entropy cocycle of (M, <p, ui) is c* = — log |det D(j) t \ | M and 
the entropy production observable is given by 



u = - - log | detity*! 



= -divX| 

M,t=0 



M ' 



Assumptions |(E 1 )|(E2)| are clearly satisfied. 

A special class of differentiable dynamics is provided by the so called Gaussian thermostated systems. Consider a 
Hamiltonian system with n degrees of freedom. The phase space of the system is R™ © R™ (or more generally the 
cotangent bundle of a smooth manifold). For simplicity, we shall assume that its Hamiltonian H is C 2 and that the 
finite energy subsets {(p, q) | H(q,p) < E} are compact. These assumptions ensure that the equations of motion 

pt = -VgiTOt, qt) , qt = V p H(pt, q t ), 

define a global C 1 Hamiltonian flow cftjjip, q) — (pt, qt) which preserves the energy, H o <fi H = H, and Lebesgue 
measure on R™ © R n (Liouville's theorem). 

To drive this system out of equilibrium, an external non-Hamiltonian force F(q) is applied. To prevent it from 
heating up, the energy supplied by this force is removed by a thermostat (the so-called Gaussian thermostat). This 
leads to the modified equations of motion 

Pt = -V q H(pt, q t ) + F(q t ) - e(q t ,p t ), qt = V p H(p t , q t ), 
where the thermostating force is given by 

F(q)-V p H(p,q) v7 
- Q(P ' g) = - \V p H(p,qW V ^ q) - 

One easily checks that the flow cj> generated by this system satisfies H o <f}* = H and therefore preserves the total 
energy. This flow, however, does not preserve Lebesgue measure and the entropy production observable 

a = V p • Q(q,p), 
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measures the local rate of phase space contraction. 

Fixing E € Ran H we see that Gaussian thermostated systems are special cases of differentiable dynamics with 

M = {{p, q) G K™ © R n | H(p, q) = E}. 

Two other well-known thermostating mechanisms are the isokinetic and Nose-Hoover thermostats. Models using 
these thermostats have been constructed to describe various phenomena like shear flows |ECM, CL|, heat conduc- 
tion IIHHPllPHl . and turbulent fluids IIGa31 lGRS I. They all have in common that the dynamics is described by a 
deterministic finite-dimensional dynamical system on a compact manifold and are very convenient for numerical 
studies llHollEMllDol . 

A well-known model in this class is a Sinai billiard with an external electric field MCELS 1 1 ICELS21 IChll ICh2l 
lYollRYL General mathematical results concerning thermostated Hamiltonian models can be found in EGCl||GC2l 
lGaUlGi2llRu2llRu3URu4llRu5l 



3.6.1 A micro-canonical ideal gas out of equilibrium. 

In this section we consider an exactly solvable thermostated system - a gas of A > 1 identical, non-interacting 
particles moving on a circle. The phase space is R N x T N equipped with Lebesgue measure and the Hamiltonian 
is 

H{L,9)= l -\L\\ 



The flow <fi is generated by the system 

L 3 =F-\(L)L 3 , 



(j = l,...,N). (3.20) 



Here, F € R denotes the constant strength of the external forcing and —X(L)Lj is the thermostating force, 

N -, N 



KL)=F i -, i=^L k , u=^J2 L l- (3-21) 

k=l k=l 

The mean kinetic energy per particle u is constant under the flow <f> and we consider the dynamical system 
(M, </>, w) where 

M = {(L,0)\u = e}~S N - 1 xT N . : 
for some e > and lu is the normalized micro-canonical measure 

N 



"(f) = i[ f(L,e)S(u-e)f[dL j de j 
The Cauchy-Schwarz inequality ensures that £ 2 < u on M and the observable 
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is well defined. One derives £ = \i where fx = Fe is a constant. It follows that I — -y/eth.£ which, once 
inserted in d3.2U , allows to integrate the equations of motion d3.20t to obtain 

Ojt = 9jo + ~B-{ L jO ch^o - Vesh£ )(arctg(sh£ t ) - arctg(sh£ )) + ^ log (3.22) 
b b cn£ t 

ch£ t 

where £ t = £ + fit. The entropy production observable is 

u= (N- l)-=th£. (3.24) 

Assumptions |(El)|(E2)| are satisfied in this model. Moreover, the map d(L 1 9) = (-L, 9) is a time reversal for the 
flow and the measure w is lTRII 

One easily computes 

Since the distribution of £ induced by the measure u is 

^o)- 7a ^ 3 j£/(o(* { )-<*-'>* 

we conclude that 

«<°> - 108 (vans £(= h £)-<»-"<-"'(c..« + «) . (3^25, 

The validity of the finite time ES-symmetry, ensured by Proposition l3.5l can be explicitly checked by noticing that 

/oo 
(ch(£ - M t))-( JV - 1 )( 1 - a )(chO" (JV - 1)a d£ 
-oo 
oo 

(ch(-£ - ^))-( jv - i )( i - Q )(ch(-e))-( jv - i )« de 

-oo 
oo 

(ch(£ + ^))-( Ar - 1 )( 1 -«)(ch£)" (A '" 1)a d£. 

J —oo 

This example continues in Sections [5.5.1l and l6.4. II 

4 Thermodynamics 
4.1 Basic notions 

Suppose that our dynamical system (M, cf>x, ujx) depends on some control parameters X = (Xi,--- ,Xn) 6 1". 
One can think of the Xj's as mechanical or thermodynamical forces (affinities in the language of non-equilibrium 
thermodynamics) acting on the system. When dealing with such families of systems, we shall always assume that 
|(F1 )|(F2)] |(C)| and |(E 1 )| hold for each system (M, <f>x>&x)- The entropy production observable of (M, <f>xi^x) 
is denoted ax- We shall also assume: 
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(Tl) wo is </> invariant. 

We will write = <fi Q , lu = us , ui M = ui t , etc, and refer to the value X — as equilibrium. Under assumption 
|(Tl)| the entropy cocycle satisfies c* = c = for all t G K and consequently a = ct = 0. 

Definition 4.1 We call a family of vector-valued real observables 4? x = (^X'"" '^X^)> ^ a fl ux 
relation if for all X, 

N 

a x =X-^ x =Y J X^% ) . (4.26) 

i=i 

In what follows our discussion of thermodynamics concerns a family of quadruples (M, <fix , &Xi &x), where $x 
is a given flux relation. In concrete models arising in physics, physical requirements typically select a unique flux 
relation <t>x (see Section l4~4l for an example). We will refer to & x as the flux (or current) observable associated 
to the force Xj . Since cto = 0, if the map X H> ax is smooth we can always pick the fluxes as 

*X = / V(Ty|y =tl xdw. 

Jo 

Remark. To simplify the notation, unless otherwise stated we shall always assume that (M, <fix,<j->x,&x) is 
defined for all X g 1^. In concrete situations (e.g. like in the class of examples introduced in Section [4~4l , 
the systems may only be defined on a restricted range of the physical parameters Xi, . . . , Xjy. This causes no 
difficulties-one can either trivially extend the range of parameters to all of R N or indicate in each statement the 
range of parameters to which they apply. 

Our second general assumption concerns time reversal. 

(T2) The dynamical systems (M, 4>x,ux) are time-reversal invariant and 

&XOl)x = -*x 
This assumption implies that = for all X. 

4.2 Finite time Generalized Evans-Searles symmetry 

(i) 

& x ° 4>x- The entropy cocycle can be written as 
c x = tX ■ S x . 

Let P x be the law of i.e., the Borel probability measure on M. N such that P x (f) — W x(/(S^)) for any 
/ E B(R N ). Let r : R N -> R N be the reflection t(s) = -s and P x = P x o r. 



Let 



^x 



(i) 

where = *x ° <Px> ®Xs 
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Proposition 4.2 IfAssumptions \(TT)\(T2)\ hold, then for any t £ R the measures P x and P x are equivalent and 




The proof of this proposition is very similar to the proof of Proposition 13.31 and we will omit it. 

We shall call the universal relation j4.27j the finite time generalized Evans-Searles (GES) identity. As for the finite 
time ES identity, one can reformulate j4.27j in terms of the Laplace transform of P x . To this end, consider the 
functional 



One easily sees that it inherits many properties of the Renyi entropy et(a). For fixed X it is a convex function of 
Y e R N which satisfies g t (X, 0) = g f (X, X) = 0. The lower bounds 



hold, and in particular gt(X, Y) > — oo. Most importantly, Proposition 14.21 is equivalent to the finite time GES- 
symmetry expressed by the following proposition. 

Proposition 4.3 Under Assumptions | (Tl )|(T2) | one has 



for any X, Y G R and any tel. 

We again omit the proof which follows the same lines as the proof of Proposition l3.6l 

For later applications we recall the following elementary result (explicit in HLS2H and implicit in MGallO which we 
shall call the symmetry lemma. We say that a function a(X, Y) is C 1,2 in an open set O C M. N x M. N if all the 
partial derivatives c^a, 9y 4 a, dy i dy i a, 9x,0y;& an d dyfixjO' exist and are continuous in O. 

Lemma 4.4 Let the function a(X, Y) be C 1,2 in a neighborhood of '(0, 0) € R^ x and such that 




(4.28) 




g t (X 7 Y)=g t (X,X-Y), 



(4.29) 



a(X,Y) 



a(X,X -Y). 



Then 



d Xh dy ia (X,Y)\ x=Y=0 



-d Y] d Yk a(XX)\x=Y= . 



Proof. The identity 



d Xk a(X, Y) | X=Q = d Xk a{X,X-Y)\ x =a = (d Xk a) (0, -Y) + (dy, a) (0, -Y) 



leads to 



d Yj d Xk a(X,Y)\ x=Y=0 = -d Yj d Xk a(X,Y)\ x=Y=0 - d Yj d Yk a(X,Y)\ x=Y=0 . 
The equality of mixed partial derivatives d Yj d Xk a — d Xk d Yj a implies the statement. □ 
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4.3 Finite time linear response theory 

For any real or vector valued observable / we set 

(f)t = ^J w(f.)d*. 

Finite time linear response theory is concerned with the first order perturbation theory w.r.t. X of (3»x)t- Hence, 
in addition to |(Tl)f(T2)| we assume: 

(T3) The function X i-> ($ x )t is differentiable at X = for all t. 

The finite time kinetic transport coefficients are defined by 

Ljkt = d x h ($x)t\x=a- 

Since 

(a x )t = X ■ (<f> x )t = X o x kLjkt + o{\X\ 2 ) > 0, (4.30) 

the real quadratic form determined by [Ljkt] is positive semi-definite. This fact does not depend on the TRI 
assumption |(T2)| and does not imply that Ljkt — Lkjt- We shall call the relations 

Ljkt — Lkjt, 

the finite time Onsager reciprocity relations (ORR). As general structural relations, they can hold only for ITRII 
systems. 

Another immediate consequence of Equ. ( 14.301 ) is: 

Proposition 4.5 Suppose that \(Tl)\ holds and let <t»x, &x be two flux relations satisfying |(T3)| Then the corre- 
sponding finite time transport coefficients satisfy 

Ljkt + Lkjt = Ljkt + Lkjt- 

If the finite time ORR hold, then Ljkt — Ljkt- 

In the next proposition we shall show that the finite time ORR follow from the finite time GES-symmetry estab- 
lishing along the way the finite time Green-Kubo formula. 

Proposition 4.6 Suppose that \(TY%(T2)\ hold and that the function gt(X, Y) is C 1,2 in a neighborhood of (0,0). 
Then \(T3j\ holds and: 

(1) The finite time Green-Kubo formula holds, 




(4.31) 



(2) The finite time Onsager reciprocity relations hold, 

Ljkt = Lkjt- (4.32) 
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Proof. From the definition ( 14.28b we derive dy j gt{X, Y) | y _ Q = — £($^) t , and hence 

= ajc*<*^>t| x=0 = -^x fe ^,5t(x,y)| x=Y=0 . 

The finite time GES-symmetry and the symmetry lemma yield 

L m = j t d Y] d Yk g t {XiY)\ x=Y=0 = 1 jf J u>($W<f>U) )dsids2 . 

Assertion (2) follows from the equality of mixed derivatives dy } dy k gt — dy k dy^ gt- Since to is invariant, we further 
get 

which proves Assertion (!.).□ 

We finish this section with two remarks. 

Remark 1. The identity ($( fe )$i j) ) o $ = $( fe )$f?] implies that 

L ifet = ^*o;($W$W)(l-|)d S . 

Remark 2. The covariance matrix D t = [Djfct] of the vector- valued random variable 

1 I* 

Til * sds ' 

with respect to ui is 




(4.33) 



The time-reversal plays no role in ( 14.331 1. However, if the assumptions of Proposition 14.61 hold and L t = [Ljkt], 
then obviously 

D t = 2L t . 

These are the finite time Einstein relations which link the finite time covariance of fluxes in equilibrium to the 
finite time kinetic transport coefficients. Together with Proposition l4.6l thev constitute the finite time Fluctuation- 
Dissipation Theorem. We shall return to this topic at the end of Section \5?2\ 

4.4 Example: Thermally driven open systems 

We consider a system S, with phase space Ms = K™ s © K ns and Hamiltonian Hs(ps, Qs), coupled to N heat 
reservoirs Ri, . . . , R^. The phase space and the Hamiltonian of the j-th reservoir are Mj = M. nj © M. nj and 
Hj(pj,qj). The phase space and the Hamiltonian of the composite system are 

M = M s © Mt © • • • © M N , H Q (p, q) = H s (ps, qs) + fTi(pi, qi) + • • • + H N (p N ,q N ), 



and we denote by m the Lebesgue measure on M. 
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The coupling between the system S and the j-th reservoir is described by the Hamiltonian Vj(ps,Pj,qs, Qj)- The 
full Hamiltonian is 

N 

H(p,q) = H (p,q) + V(p,q) = H (p,q) + ^2v j (ps,p 3 ,qs,qj)- 

3=1 

We assume that H is C 2 and that the finite energy subsets {(p, q) | H(q,p) < E} are compact. These assumptions 
ensure that H generates a global Hamiltonian flow of class C 1 on M. For any C 1 observable F, 

where { • , • } denotes the Poisson bracket, {F, G} = V ' q G ■ V P F - V P G • V q F. 

The state of the combined system in which each reservoir is at thermal equilibrium at inverse temperature f3j and 
the system S at inverse temperature j3 is the product measure 

I c -^ s -£f =1 m 

Introducing the control parameters Xj ■ — (3 — /3j, we can rewrite it as 

lux = ^e~ pHo+?: ?=i XjHi m. (4.34) 
Z 

The dynamics does not depend on X and we set (f> x = </>*. Note that u is not invariant under the flow 0*. In order 
to satisfy hypothesis |(T 1 )| one modifies (14.34b as 

lu x = -c-^ Hs+v ^U M m = le-^+E^i W m . (4 .35) 
Z Z 

With this definition, loq is the Gibbs canonical ensemble at inverse temperature j3 and is invariant under <p l . More- 
over, if the reservoirs have a large spatial extension and the coupling Hamiltonians Vj are well localized, the states 
( 14.34b and ( 14.35b describe the same thermodynamics. 

For the reference state ( 14.351) the entropy cocycle is given by 

N N ,t 



= *i (Hit -H i ) = -Y i X i J ^H js ds, 



i'=l 3=1 

and we have the flux relation 

N N 
3=1 3=1 

The flux observables do not depend on X and are given by $0') = {Hj,V}. Since 

Hjt -Hi = - f ^ds, 
Jo 

the flux observable $W describes the flow of energy out of the j-th reservoir. The time reversal in physical systems 
is usually given by the map d(p, q) — (— p, q) and the system (M, <f>x,^x) is then lTRII provided H o d = H. 

We shall investigate a simple example of thermally driven open system in the remaining part of this section, to be 
continued in Sections [5. 5. 21 and |6. 4. 21 
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4.4.1 The harmonic chain 



Quadratic Hamiltonians provide instructive examples of open systems whose non-equilibrium characteristics can 
be computed in a closed from |LS1 1. From the mathematical point of view they are special cases of the Gaussian 
dynamical systems discussed in Section [9] Since the entropic fluctuations of such models are studied in detail in 
the forthcoming paper [JLTP|, for reasons of space we shall be brief. 

For a finite subset T = {n, n + 1, . . . , n + k — 1} C Z we set M T = R r © R r and we define 

M r ^ M h> H*<p, q ) = £ ^ + (fe " 2 fe - l)2 , 

where we set p x — q x — for x $ T (Dirichlet boundary conditions). 

For some integer m > 0, let Ts = {— m, . . . , m} and set Hs = H Ts . This Hamiltonian describes a finite 
harmonic chain. We shall couple it to two large heat reservoirs, Rl and Rr, at its two ends. For this purpose, let 
d> m and set Tl = {—n, . . . , —m — 1}, = {m + 1, . . . , n}. The Hamiltonians of the two reservoirs are 

H L =H r \ Hr = H Fn . 



The Hamiltonian of the composite (but still decoupled) system is 

H = H Tl + H Ts + H Vr =H l + H s + H r . 
Finally, define the Hamiltonian of the coupled system by 

H = H 1 



The coupling Hamiltonian is given by 

V = H - H = Vl + Vr = Q — m — lQ—m ~ imlm+U 

and is independent of n. Since the equations of motion induced by Hq and H are linear, the associated Hamiltonian 
flows are linear group which we write as e tCo and e tc respectively. 

Let us denote by h, h^, Iir the real symmetric matrices corresponding to the quadratic forms 2H, 2Hl, 2Hr. The 
reference state ujx is the centered Gaussian measure of covariance 

D x = ((3h-k(X))~ 1 , 

where 

k(X) = X L h L X R h R . 

For /3 > the set Op = {X £ M. 2 \ /3h — k(X) > 0} is an open neighborhood of 0. The dynamical system thus 
obtained is well defined for X £ Op, is lTRIl and clearly satisfies Assumptions |(T1 )| and |(T2 )| 

It is a simple exercise in Gaussian integrals to show that Assumption |(El)| is satisfied for X small enough. However, 
note that the flux observables 

$W = {H L ,V} = -p-m-lQ-m, 
^ = {Hr,V} = - Pm +iq m , 

as well as entropy production ax = Xl& l ^ + Xr<&( r > are unbounded. Thus, Hypothesis |(E2)| is not satisfied. 
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Propositions 14.21 and 14.31 apply. Moreover, the functional g t (X, Y) reduces to a Gaussian integral that can be 
computed explicitly 

g t (X,Y) = -ilogdet (/- D x (e tc 'k(Y)e tc - , 

with the convention that log.T = — oo for x < 0. Since the groups e tc is uniformly bounded, the function 
(X, Y) i — ^ gt{X, Y) is real analytic on an open neighborhood of (0, 0) in Op x R 2 which is independent of t G K 
and Proposition ^. 61 applies too. 



The validity of the finite time GES-symmetry, ensured by Proposition 14. 31 can be explicitly checked as follows. 
Energy conservation, e tc he tc — h, yields 

D x l - (e tc *k{Y)e tc - k{Y)) =/3h- k{X — Y) — e t£ * k(Y)e tc 

= e tc ' (j3h - k{Y) - e~ tc ' k{X - Y)e~ tc ^ e tc 

= e tc * {Dx 1 - {^ tC 'k{X - Y)c- tc - k(X - F))) e t£ . 

This relation and Liouville's theorem, dete tc = dete t£ = 1, imply that g t (X,Y) — g^ t (X,X — Y). Finally, 
UUyields g- t {X, X - Y) = g t {X, X — Y) and the finite time GES-symmetry follows. 



5 The large time limit 
5.1 Entropy production 

For any observable / we set 

(/> + 

whenever this limit exists. 

In this section in addition to |(El)[|(E2)| we assume: 

(E3) The limit (er) + exists and is finite. 
The entropy balance equation yields the basic result: 
Proposition 5.1 (a) + > 0. 
We shall say that the dynamical system (M, (j>, u>) is entropy producing if it satisfies 

(EP) > 0. 

The validity of |(E3)| and |(EP)| are dynamical problems that can only be answered in the context of concrete models. 
In this section we shall discuss several structural results which shed some light on these central issues. 

Proposition 5.2 Suppose that 

(a) t = (a)+ + 0(t- 1 ), (5.36) 
as t — > oo. Then (cr)+ = implies that there exists v € Si (1 Af u satisfying Ent(^|cij) > — oo. 



= lim (f) t , 

t— too 
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Proof: We shall use the properties of relative entropy listed in Theorem |2.1| Suppose that (a) + = 0. The entropy 
balance equation and ( 15.361 1 yield 



|Ent(wt|w)| 



(w(ct s ) - (<r)+)ds 



o 



O(l) 



Hence, there is C such that Ent(wt|w) > C for all t > 0. Set £>{(•) = J CJ s (-)ds. The concavity and the 
upper-semicontinuity of the relative entropy yield Ent(d) t |w) > C. By compactness, there exists v G S and a net 
t a — > oo such that — > v. It follows easily that v 6 <S/ and the upper-semicontinuity implies Ent(z/|w) > — oo. 
□ 

Proposition 5.3 Let i/£5;fl Af u . Then i/(a) = 0. 



Before proving Proposition l5.3l we need a preliminary result which is of independent interest. In what follows we 
shall say that a sequence t n t oo is regular if J2 n e~ at ™ < oo for all a > 0. 

Lemma 5.4 Let t n be a regular sequence. Then, for ui-a.e. x, 

1 /"*" 1 /**» 

liminf — / a s (x)ds > 0, limsup — / er_ s (a;)ds < 0. (5.37) 

n— >oo t n Jq n — >oo t n Jq 

Proof. We will prove the first relation in (15.37b . a similar argument yields the second. 

Let X n = t^ 1 J*" a s ds and A = {x € M | liminf n _>oo X„(a;) < 0}. We need to show that lj(A) = 0. Since 
A = Uk>iAk with Ak = {x € M | liminf n^oo X n (a;) < — 1/fc}, it suffices to show that w(Afc) = for all 
integers k > 1. Set A* = A Ut | u and note that cj(A _t ) = w_ t (l) = 1, for all t. The Markov inequality gives 

uj({x e M| A~*(a;) > A}) < A -1 , 

for A > 0. Since A _t " = e~*™ Xn , we have 

e M | JC„(a:) < -a}) = cj({x G M | A _ *»(a:) > e Qt "}) < e" *". 

Hence, 

€ M| JC n (a:) < 1/A;}) < oo, 

n 

and the Borel-Cantelli lemma yields that ui(Ak) = 0. □ 

Proof of Proposition |53J By the ergodic theorem, there is a £ L 1 {M, dv) such that 

i r 

lim — / a± s (x)ds — a(x), 
n->oo n J 



for i^-a.e. x. Since v is invariant one has v(cr) = v(cr). Since v is normal w.r.t. oj, Lemma I3~4l implies that 
v(o) = and the statement follows. □ 

Corollary 5.5 Suppose that Si n Af u =/= and that there exists a sequence t n f oo smc/i f/zaf 



for uj-a.e. x. Then (<x) + = 0. 



lim — / a s (x)ds = (<t) + , 
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Proof. Let v e <Sj n M u . Then 



1 '*« 



(cr) + = j/ lim — / <7 s (x)ds = f(<x) = 0. 



□ 



The results of this section establish that, under very general conditions, the dynamical system (M, 4>, uS) is entropy 
producing iff Si n Af u = 0- 

5.2 Linear response theory 

Consider a family (Af, (j>x,^x, &x) satisfying |(T 1 )| \(TTj\ and [03)1 In this section we are interested in the large 
time limit of the time averaged expectation values of individual fluxes and validity of the linear response theory. 
In addition to (T1)-(T3) we assume: 

(T4) The limit (<&x}+ exists for X small enough and is differentiable at X = 0. 

(T5) u}(^^ j) ) = 0{t~ l ) as t ->■ oo. 

The verification of (T4) and (T5) is a dynamical problem that can be answered only in the context of concrete 
models. |(T2)1 and |(T5)1 imply that = 0(* -1 ) as i -> -oo. 

The kinetic transport coefficients are defined by 



L j k = d Xh {*% ) )+\x=o' < 5 - 38 ) 



Since (<Jx)+ = Y^j Xj ($ < x > )+ — ^> tne rea ^ quadratic form determined by [Ljk] is positive semi-definite. 
The kinetic transport coefficients satisfy the Onsager reciprocity relations (ORR) if 

L jk =L kj , (5.39) 

and the Green-Kubo formula holds if 

L jk = - / W (4>( fc )$«)ds, (5.40) 
where, unless otherwise specified, = lim^^ /_ ( . Note that d540l 0391 

The finite time linear response theory leads to a natural axiomatic program for the verification of (15.39b and ( 15.401 ) 
based on the following 

Proposition 5.6 Suppose that (T1)-(T5) hold. Then the following statements are equivalent: 



(1) The Green-Kubo formulas ( 15.401 ) hold. 

(2) lim t _j. 00 Ljkt = Ljk- 
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Proof. Set 



i ft 

F(t) = - / w{& k) $i j) )ds, 



and notice that 

L jkt = 1 1* W ($W$C7)) (i_ ds = \j\{s)As. 

By the fundamental property of Cesaro's mean, (1) =>• lim^oo F(t) = Ljk => (2). On the other hand, Hypothesis 
|(T5)| and Hardy-Littlewood's Tauberian theorem (see e.g. [Ko|) yield (2) =>■ lim^oo F(t) — Ljk (1). □ 



We would like to add several remarks regarding Proposition^ 

Remark 1. Finite time linear response theory requires the minimal regularity assumptions (T1)|(T3)1 and in 
particular no ergodicity assumption. It is valid in practically all models of interest. Assumption |(T4) states that the 
basic objects of linear response theory are well-defined (existence of the Ljk 's) and is of course necessary to have 
a meaningful theory. Condition (2) of Proposition |52] can be reformulated as 



8 Xk ( Urn (*g>) t ) = lim (d Xk ) , (5.41) 

Vt->oo / x=0 t—Hx \ X=0/ 

i.e., as an exchange of the two limits t — > oo and X — > 0. Even though the existence of the improper integral 
in d5.40t does not require any decay of the correlation function t t-> w($^ fe - ) $^), Assumption |(T5)| provides the 
minimal decay assumption which ensures that 



1 r* 

lim L jkt = lim - / w($ (fc) $W)ds 

;— >cx) t— voo 2 J 



t— f oo 

Note however that Assumption |(T5)| is needed only for the Green-Kubo formula and that Condition (2) auto- 
matically implies ORR. Assumptions |(T4)| and |(T5)| are ergodic in nature and are typically difficult to verify in 
physically interesting models. A particularly delicate aspect is differentiability of the function X M> (3>x) + - 

Remark 2. The proposed program for the derivation of linear response theory is in a certain sense minimal. On 
physical grounds one would like to have an additional estimate 

(*x)t = LX + Er(X,t), (5.42) 

where the error term satisfies 

x ~>0 t>0 \x\ 

with the rate of convergence/range of parameters that allow to draw physical/numerical conclusion from d5.42t . 
This point is related to van Kampen's objections against linear response theory [Ka KTH|, see |CELS2| for a 
discussion. 

Remark 3. In some models the following well-known result (the multivariable Vitali theorem) can be effectively 
used to verify (EH) (see HEQ). Let I e = {X e 1* | \X\ < e} and D e = {X e C N | \X\ < e}. 

Proposition 5.7 For allt > let Ft ■ D e — > C be an analytic function such that 

sup \F t (X)\<oo, 
xeD c ,t>o 

and assume that 

lim F t (X) = F(X), (5.43) 

t— foo 

exists for X € I e . Then the limit ( 15. 43 i exist for all X £ D e and is an analytic function on D e . Moreover, as 
t — > oo, all derivatives of F t converge uniformly on compact subsets of D e to the corresponding derivatives of F. 
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We shall point out some mathematical intricacies regarding the interchange of the limit and derivative in ( 15.411 ) on 
a simple example in Section [8731 



Our final topic in this section is the Fluctuation-Dissipation Theorem (recall its finite time counterpart discussed in 
Remark 2 of Sectionl4~3l. 



Definition 5.8 Suppose that (T1)-(T4) hold. We shall say that the Fluctuation-Dissipation Theorem holds for 
(M, </>x,w x ,*x)) if: 

(1) The Green-Kubo formulas (and hence the Onsager reciprocity relations) hold for the kinetic transport coeffi- 
cients L = [Ljk], 

(2) The Central Limit Theorem holds for <& = ($ W : . . . j (j>W ) w.r.t. [M, <fi, oj) with covariance matrix 

D = 2L, (5.44) 

i.e., for any Borel set B C M. N , 



lim oj ( <j x e M 

t— > oo 



where /id is the centered Gaussian measure of covariance D on R". 
Remark. The celebrated Einstein's relations ( 15.44b link equilibrium fluctuations to kinetic transport coefficients. 



Just like Proposition l5.61 the Fluctuation-Dissipation Theorem is "forced" by its universally valid finite time coun- 
terpart. With regard to the proof of the Central Limit Theorem, we mention the following result of Bryc (Bryj. D e 
and I e are as in Proposition l5.7l 

Proposition 5.9 Assume that \(Tl)\ holds. Suppose that for some e > the function 

g t (0,Y)=logu>(e- Y -ti** ds ) ) 

is analytic in D e , satisfies 

sup -\g t (0,Y)\ < oo, 

Y£D € ,t>l t 

and that 

lim i fft (0,y), (5.45) 

t— >oo t 

exists for all Y € l t . Then the Central Limit Theorem holds for $ = ("I^ 1 ), • • ■ , w.r.t. (M, 4>,uj) with 

covariance matrix 

D jk = t lim j w(*W#W)) ^1 - ds. (5.46) 



Remark. If ^> g L°°(M, du>), then the function Y gt{0,Y) is real analytic. The location of the complex 
zeros of the entire analytic function Y M> u;(e -r '-'° * sds ) determines the region of complex plane to which 
<?t(0, Y) extends analytically. If in addition ["(T2)||(T3)| and |(T5)| hold, the existence of the limit ( 15.461 ) implies that 

i r°° i 

ljmL jkt = - / u @m*U))ds = -D jh . 
t^oo 2 J_ co 2 

The Fluctuation-Dissipation Theorem is the pillar of non-equilibrium statistical mechanics in the regime where the 
thermodynamic forces are weak. The far from equilibrium case is discussed in the next section. 
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5.3 The Evans-Searles Fluctuation Theorem 



We start by recalling some basic facts of the Large Deviation Theory (see e.g. MDZlER ). 

Definition 5.10 A vector-valued observable f = (f^\ ■ • • , /^) satisfies a Large Deviation Principle w.r.t. to 
(M, 4>,uj) if there exists an upper-semicontinuous function I : M. N — > [—ex, 0] with compact level sets such that 
for all Borel sets G C R N we have 

sup I(Z) < lim inf - log cj 

zed 1 

< lim sup - log u> 
where G denotes the interior of G and G its closure. 



x e M 

x £ M 



i jT f s (x)ds£G^ 

- / f s (a:)ds eC?}) < sup/(Z). 
t Jo ) J zee 



(5.47) 



The following standard result goes under the name of Gartner-Ellis Theorem and will be used repeatedly, see e.g. 
llEIllDZl for a proof. 

Proposition 5.11 Assume that the limit 

h{Y)= lim -bgwfe-Jif^W-) 
exists in [— oo, +oo] for all Y £ M. N and is finite for Y in some open neighborhood of £ R*. 

(1) Suppose that h(Y) is differentiate atY — 0. Then, the limit 

(f)+= lim ~ fujiQds, 
exists and (f) + = — V/i(0). Moreover, for any regular sequence t n one has 

lim 1 / " f s (x)ds = <f)+, 

™^°° tn Jo 

for ui-a.e. x. 

(2) Suppose that h(Y) is a lower semicontinuous function on M. N which is differentiable on the interior of the set 
V = {Y £ R N | h(Y) < oo} and satisfies 

lim |Vft(y)| = oo, 

for all Yo € dT>. Then the Large Deviation Principle holds for f w.r.t. (M, (f>, lu) with the rate function 

I{Z) = inf (Y-Z + h(Y)), 

YeR N 

i.e., —I(Z) is the Legendre transform ofh{—Y). In particular, I{Z) is concave. 

Remark 5.12 The conclusion of Part (2) holds in particular if h(Y) is differentiable on K . There are other 
(local) versions of Gartner-Ellis theorems that are useful in applications. Suppose, for example, that the function 
h(Y) is finite, strictly convex and continuously differentiable in some open neighborhood B C of the origin. 
Then Part (1) holds and a weaker version of part (2) also holds: the large deviation bounds ( 15.47b hold provided the 
set G is contained in a sufficiently small neighborhood of the mean (f) + (see Lemma XIII. 2 of HHH1 and Section 
4.5 of||DZ)). 
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Let (M, 0, w) be a lTRII svstem. Recall that 

e t (a) =logw(e- Q -/"o ff . d «) . 

We suppose: 

(ES) The Evans-Searles functional (ES -functional for short) 

R3a4 e(a) = lim -et(a) e [—00,00], 
t— ►oo £ 

exists. 

Propositions 13 .4l and l3 .5 1 yield the basic properties of the ES-functional : 

Proposition 5.13 (1) e(a) is a convex function of a. 
(2) It satisfies the ES-symmetry 



(3) 



e(a) = e(l - a). (5.48) 

e(0) = e(l)=0, 

e(a) < jfa6[0,l], 
e(a) > otherwise. 



(4) /f satisfies the lower bound 



with 



1 

a-- 



e(a) > 

£+ = limsupw(E*). 



S+ (5.49) 



/n particular, i f\(E3)\ holds then E + = (cr) + . 

We emphasize that the ES-symmetry ( 15.48b is an immediate consequence of the finite time ES-symmetry. 
Using Proposition l5.1 ll we obtain 

Proposition 5.14 (1) Suppose that e(a) is differentiable at a — 0. Then \(ETj\ holds and (o~)+ = — e'(0). 

(2) Suppose that e(a) is differentiable for all a £ K. Then the Large Deviation Principle hold for the entropy 
production observable a w.r.t. (M, <p, uj) with the concave rate function I(s) — inf a6 R(sa + e(a)). Moreover 

I(s)=s + I(-s). (5.50) 

Proof. We only need to prove d5.50t . Using d5.48t we have 

I(s) = inf(sa + e(a)) = inf(s(l — a) + e(l — a)) — s + inf (— sot + e(a)) = s + I{— s). 

a a a 

□ 

The relation d5.50t is called the ES-symmetry for the rate function I(s). 

Consider now a family (M, <px,^x, &x) indexed by X 6 R N and satisfying Assumptions |(T1 )| |(T2)| and |(T3)| 
We assume: 
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(GES) The Generalized Evans-Searles functional (GES -functional) 

g{X,Y)=Km ]gt(X,Y), 

exists for all X, Y. 

g(X, Y) is a convex function of Y and the finite time GES-symmetry implies that 

g(X,Y)=g(X,X-Y). 
We shall refer to this relation as the GES-symmetry. 

Proposition 5.15 (1) Suppose that Y H> g(X, Y) is differentiable at 0. Then \(T4)\ holds, 

(<t> x ) + = -V Y g{X,Y)\ Y ^ 

and ift n is a regular sequence, then 

1 /•*" 

lim — / $ Xs (x)ds = (*x)+> 

n-i-oo t n J 

for LOx-a-e. x. 

(2) Suppose that Y g(X, Y) is differentiable for all Y. Then the Large Deviation Principle holds for the flux 
observables w.r.t. (M, (fix ,Wx) with the concave rate function 1 x (s) — m fYeR N (^' s + .9(^?^))- Moreover 

I x (s) = X-s + I x (-s). 

(3) Suppose that g(X, Y) is C 1,2 in a neighborhood of (0, 0). Then the kinetic transport coefficients are defined 
and satisfy the Onsager reciprocity relations. 

(4) In addition to the assumption of (3) suppose that \(T5)\ holds and that for some e > 0, 

sup i| 5t (0,y)|<oo. (5.51) 

YeD e ,t>i t 

Then the Fluctuation-Dissipation Theorem holds for (M, (fix, u) x , 

Proof. (1) and (2) are immediate from Proposition 15 . 1 1 1 and the symmetry of Ix(s) is proved as in Proposition 
EH 

(3) By Assertion (1), ($^) + = dy^giX, Y) \ Y _ Q , hence the GES-symmetry and the symmetry lemma yield 

L jk = dx k d Y] g(X,Y)\ x=Y=Q - -\d Yi d Yk g(X,Y)\ x=Y=0 . (5.52) 

Since the partial derivatives on the right hand side are symmetric in j, k, we have Ljk — Lkj. 

(4) From (I5.52l > and the fact that 

d Y] d Yk \\ogg t (Q,Y)\ Y=0 = - (l- y) d», 
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we see that the Green-Kubo formula holds iff the limit and the derivative in the expression 

d Y] d Yk g(0, Y)\ Y=0 = d Y] d Yk Km ~ 9t (0, Y)\ Y=Q , (5.53) 

can be interchanged. This is ensured by the assumption (15.511 ) and Proposition 15.71 Similarly, Proposition 15.91 
yields the CLT. □ 

We shall say that a given ITRll model satisfies the Evans-Searles Fluctuation Theorem if the respective functionals 
e(a)/g(X,Y) exist and are differentiable/C 1,2 . It follows from Propositions 15 . 1 41 and 15 . 1 5l that the the Evans- 
Searles Fluctuation Theorem can be interpreted as an extension of the fluctuation-dissipation to the far from equi- 
librium region. 

Remark 1. The sufficient and necessary condition for the validity of the Green-Kubo formula is that the limit and 
the derivative in the formula ( 15.531 ) can be interchanged. Assumption ( 15.511 ) provides a convenient criterion for 
validity of this exchange which will be satisfied in several examples that we will consider. In general, however, 
there may exist other mechanisms that will lead to the justification of ( 15.531 ), see for example the proof of the 
Fluctuation-Dissipation theorem for the Sinai billiard with small external forces in [ CELS 1 CELS2 Chll ICh21 . 

Remark 2. It is instructive to compare (3) and (4) to the finite time based derivation of the linear response theory 
presented in Section l5T2l 

Remark 3. In some models where the entropy production observable is unbounded the ES-functional e t (a) is 
finite only on an open interval containing [0, 1]. In this case one can still formulate a meaningful Evans-Searles 
Fluctuation Theorem, see Section[9]for an example. The same remark applies to the GES -functional (see Section 
13321) . 

5.4 The resonance interpretation of e(a) 

Under suitable regularity conditions, the identity 

e t (a) =log(l,e*M), (5.54) 

for p — 1/a leads to identification of e(a) with a resonance of L p -Liouvillean L p . In this section we state two 
general results regarding this identification. 

Since sp(L p ) C {z | |Rez| < m p }, the resolvent (z — L p )~ x is a well-defined operator valued function analytic 
in the half-plane Re z > m p . Note that |e(a)| < m p and that 

/>oo 

(l,(z-L p )- 1 l)= / e e ^~^- tz dt, 
Jo 

for Kez > m p . 

In the next two propositions a G lis fixed and p = 1 /a. 

Proposition 5.16 Suppose that for some 7 > and c£R, 

e t (a) =te(a)+c + 0(e- 7t ), 

as t — > 00. Then the function z > (1, (z — has a meromorphic continuation from the half-plane 

Re z > m p to the the half-plane Re z > e(a) — 7 and its only singularity there is a simple pole at z — e(a) with 
residue e c . Moreover, for any e > and j 6 {0, 1}, 

sup f \{l,{x + i y -L p y- 2 l)Y +1 dy < 00. 

x>e(a)— 7+e J \y\>e 
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This result has the following converse: 

Proposition 5.17 Suppose that the function z i— > (1, (z — L p )~ l l) has a meromorphic continuation from the half- 
plane Re z > m p to the half-plane Re z > e(a) — "/for some 7 > and that its only singularity there is a simple 
pole at z — e(a). Suppose also that for some e > and any j € {0, 1}, 

sup f \(l,(x + iy- L p y- 2 l)y +1 dy < 00. 

x>e(a)— 7 J\y\>e 

Then 

e t (a) =te(a)+c + 0(e- 7 *), 

as t — > 00. 

The proofs of Propositions l5. 16l and l5. 17l are standard (see [ JP3 1) and for reasons of space we will omit them. 
After introduction of a suitable transfer operator the resonance interpretation of g(X, Y) is very similar. 
In the discrete time case (recall d3.19t ) instead of the resolvent [z — L p ) _1 one considers 

00 

ft(z) = X>— (i,c/;i). 

71=0 

Propositions 15 . 1 61 and 15 . 1 71 hold in the discrete case after obvious modifications. 

The reader familiar with classical results in spectral theory of Ruelle transfer operators MBo2l IRull IBall IBKLI 
GL1 , GL2 Ba2 , Ba3 1 might be surprised at our insistence on the Hilbert space framework. It is however precisely 
in this framework and through the link with Tomita-Takesaki theory |BR] that Ruelle transfer operators (L p - 
Liouvilleans) naturally extend to the non-commutative setting. The Banach space framework, which is dominant 
in the classical presentations, emerges through the complex spectral deformation technique which is a natural tool 
to study resonances of the L p -Liouvilleans in the non-commutative setting (see MJP3II for the case p — 00). 



5.5 Examples 

5.5.1 The microcanonical ideal gas 

In this section we investigate the large time limit in the example of Section [3.6.1| For F 7^ it follows from Equ. 
( T3T24j ithat 

lim a t = {N-iy- 1 l, 

holds for w-a.e. (L, 9) € M, hence 

(a)+ = (TV _ i)^2 >0 . (5.55) 
The generating function ( 13.25b can be expressed in terms of the associated Legendre function as 

e t (a) = log (T(N/2) f_j ^-^-i/ai-va^*) J • 



From the asymptotic behavior for z —> +00 (see e.g. Equ. (8.766) in [GR]) 

y/TT \L (1 + n — m) 1 (—n — m) J 
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we obtain the ES-functional 



e(a) = lim ~e t {a) = -{a)+ [ ~ - 



a-- 



This function is not differentiable at a = 1/2. However, it is differentiable near a = 0, we conclude that the 
entropy production observable satisfies a (local) Large Deviation Principle with rate function 



I(s) 



1 



if Isl > (a) 



(s -(*)+) if m <<*>+, 



near its mean value (<r) + . 

With F as control parameter, we have the flux relation a = F (N — l)e -1 ' 2 th^, and hence the flux observable 

N- 1 



-thf. 



From Equ. ( 15.55b we conclude that ($) + = (N — l)e 1 / 2 sign F is not differentiable at F — and linear response 
theory fails for this model. We remark that the finite time Green-Kubo formula reads 



Of 



w($ s )ds 



F=0 



(JV - l) 2 t 
iV 2e' 



and that i t diverges as t — > ex). 



5.5.2 The harmonic chain 

We continue with the example of Section 14.4.11 We again omit the details of the calculations which the interested 
reader may find in IIJLTPI . 

Since the reservoirs in Section l4.4.1l are all finite, one has 

(® iL/R) )+ = \ J o ^x{^ s L/K) ) ds = Km Itr (d x (h L/R - e tc * h L/R e tc )) = 0. 
In particular, (a x )+ = X L (^)+ + X R (^) + =0. 

To get a non-vanishing entropy production, we must perform the thermodynamic limit of the reservoirs (i.e., take 
n — > oo keeping m fixed) before taking t — > oo. We shall not be concerned here with the existence of a limiting 
dynamical system. However, it should be clear from our discussion that the limiting dynamical system exists as a 
special instance of the Gaussian dynamical systems of Section [9] In the following, we denote the dependence on 
n of various objects of interest by the superscript . 

The phase space of the composite system has a natural embedding in the real Hilbert space T~L = ^jjj(Z) ^r(Z). 
We denote by H — Hl © Hs ®T~Lr the decomposition of this space corresponding to the partition 

Z = {x e Z | x < -m} U{x£Z \ - m < x < m}U {x £ Z\x > m}, 

and by p^, p$, p R the corresponding orthogonal projections. The operators h^ n \ h R and k^ n \X) have 
strong limits in this Hilbert space as n —> oo. We shall denote these bounded self-adjoint limits by h, h^, h R and 
k(X). For example 

s-lim/iW - h- ( 1 ° ^ 
s n ^T n n \ I- A J' 
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where A is the finite difference Laplacian on £ 2 (Z) 

(Au) x = u x -i - 2u x + u x +i- 
In the same way, the generators C , as well as their adjoints C^*, £( n )* have bounded strong limits £ , C 

s - lime* £ ° ' = e* £ °, s — lime t£< "' = e t£ , 



and Cq, £*. It follows that 



and similar relations for the adjoint groups hold uniformly on compact time intervals. 

Denote by cf> L / R the finite rank self-adjoint operator associated to the quadratic form 2<&( L / R ^ (recall that it does 
not depend on n). For Y = (Y L ,Y R ) G R 2 set 4>(Y) = Y L <f) L + Y R (j) R . It follows from 

Jo 

that 

lim e tc ^"k^(Y)e tc - k^(Y) = e tc ' k(Y)e tc - k(Y) = e s£ * <j>(Y)e sC da, 

n->oo 

holds in the trace norm for any finite tel. Since {fihW - /fc (n) (AT))" 1 is uniformly bounded and 

s - lim(/3/i (n) - fc^X))- 1 = [fih - k(X))-\ 

n— >-oo 

we conclude that 

lim {fihF* _ fcWm)" 1 fe' £(n) *fc (?1 )(y)e t£(n) - fc^F)) = - f {Bh - fc(X))- 1 e s£ >(F)e s£ ds, 

n->oo V ' Jo 

holds in trace norm. This finally yields 



ft(X,y) = lim ft w (X,y) = --logdet /+ / ((3h-k(X))- l e sL 0(F)e s£ ds , (5.56) 
with the convention that log x = — oo for x < 0. 

We are now in position to perform the t — > oo limit. The wave operators 

W±=b- \imh 1/2 e- tc e tCo h~ 1/2 (p L + p R ) = s - lim^ 1/2 e t£ * e -* £ «^ /2 (p L + p R ), (5.57) 

t— >±00 t— f±CJO 

exists and are partial isometries from Hl © Hi? to "H. The scattering matrix 5 = W+W- is unitary on Hl © 'H/i 

and commutes with the self-adjoint operator Lq = h^ 2 £oh Denoting by Hl © H R — J ffi h(A)dA the 
spectral decomposition induced by Lq and by 5(A) the fiber of the scattering matrix 5 acting on f)(A), one has 

g(X,Y)=lim±g t (X,Y) = --^ J logdet„ (A) (/-(/?/- X)- 1 (5(A)*F5(A) - ?)) dA, 

where X = XlPl + X R p R . We note that this formula remains valid for arbitrary finite harmonic system coupled 
to a finite number of infinite harmonic reservoirs, as long as the coupling v = h — ho is trace class (see |JLTP|). 
The fluxes are given by 

($«)+ = -d Yi g(X,Y)\ Y= = i- | tr h(A) ((/?/- X)- 1 ^ - 5(A)* ft -5(A))) dA, 
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which is a classical version of the Buttiker-Landauer formula (see [AJPP|). 
Explicit calculation of the scattering matrix yields the result 

fl 1 ^*,* ) = -«l0g -5 ~ — , K = 



V [fi-X L )[fi-X R ) J 1 2tt ■ 

For fixed X E R 2 such that max(Xf,, Xr) < /3, the function Y i-> F) is a real analytic in the open strip 

{y e R 2 I - 08 - x R ) <y r -y l </3- x l }. 

By Proposition l5.11l 

<*<*>>+ = -<*<*>+ ^ « (-^ - ^) ^ «(T £ - Tr), (5.58) 
where T L i R = (j3 — X L /r) ~ 1 denotes the temperature of the L/R reservoir. In particular, entropy production 

Irr \ - Y 4- Y /(b( R )\ - Z XR ^ - , S TL ~ TR ^ 2 

(*x) + -X L (<S> ) + +X R (S )+-k {i3 _ Xl){i3 _ Xr) -k TlTr • 

is strictly positive provided Tl ^ T R . By Proposition 15. Ill the flux observables (<I>( L ) , $0^)) satisfy a large 
deviation principle w.r.t. wx> with rate function 



if s L + s R ^0, 



Ix[sl,sr) 



K 



F(9) if s l = -sr = — she, 



where 

F{0) = — k 



A, = f3-(X L + X R )/2 and 8 = (X L - X R )/2. 

Writing Equ. ( 15.56b as gt{X, Y) = — ^ log det(I + At), one easily shows that the trace norm of A t is bounded by 

\\A t h<C\Y\\t\, 

while its operator norm satisfies 

\\m<c\y\, 

for any Y E C 2 and f e R with a constant C depending only on j3 and X. It follows that the bound ( 15.511 ) 
is satisfied for sufficiently small e > 0. Finally, as a consequence of the local decay estimate for the discrete 
Klein-Gordon equation 

Hypothesis |(T5)| is satisfied. It follows from assertions (3) and (4) of Proposition l5.15l that the fluctuation dissipa- 
tion theorem holds. 

The equality = — ($^') + in ( 15.581 ) is a consequence of energy conservation. More generally, for an open 

system as described in Section l4~4l energy conservation implies )+ = 0- F° r tne same reason, the rate 

function Ix (s) takes the value —00 outside of the subspace sj = and the covariance matrix D in the central 
limit theorem is singular on this subspace. One can easily avoid all these singularities by reducing the number 
of parameters. In fact, one observes that the large time characteristics of the system do not depend on the initial 
inverse temperature (3 of the small subsystem S. Indeed, the GES-functional g(X, Y) only depends on the inverse 
temperatures of the reservoirs j3 — Xj (this is a general feature of open systems). This suggests to fix the parameter 
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/? at the mean inverse temperature of the reservoirs by restricting the parameters X to the hyperplane Xj = 
of M. N . This reduces the number of parameters and consequently the number of associated fluxes by one. In our 
simple example with two reservoirs, this amounts to set Xl = —Xr, 

Harmonic systems are very special in that the study of their dynamics can be effectively reduced to an application 
of the trace class scattering theory. The more difficult case of a finite anharmonic chain coupled to infinitely 
extended harmonic reservoirs has been analyzed in |EPR1 EPR2 EH1 , EH2, IRTTllRT2ll . 



6 Non-equilibrium steady states 
6.1 Basic notions 

To discuss Non-Equilibrium Steady States (NESS) we need several additional assumptions on (M , </>, w). The first 
is: 

(NESS1) M is a complete separable metric space. 

In this case it is natural to equip S with the topology of weak convergence, i.e., the minimal topology w.r.t. which 
all the functionals 

S3u^v(f), feC(M), 

are continuous. This topology is metrizable and S is a complete separable metric space. A sequence v n € S 
converges to v iff u n (f) -> i/(f) for all / e C(M). 

With regard to Theorem l2.ll in (1) Br(M) could be replaced by Cr(M). (6) and (7) are valid as formulated except 
that in (6) the convergent nets can be replaced with convergent sequences. 

S is compact iff M is compact. More generally, a set & C S is precompact (its closure is compact) iff & is tight, 
i.e., for any e > there is a compact set K t C M such that v(K t ) > 1 — e for all v e &. 

The remaining additional assumptions are: 

(NESS2) <fi f is a group of homeomorphisms of M and the map (t, x) i— > 4> l (x) is continuous. 
(NESS3) a £ C(M). 
(NESS4) The set of states 

t 



is precompact in S. 



We denote by S+(uj) the set of limit points of &(cj) as t — > oo. S+(lu) is non-empty and uj + € <S+(w) iff there 
exists a sequence t n — > oo such that, for all / € C(M), 

1 /■*» 

lim - / u.(f)ds = w+(/). (6.59) 



n— >-oo 



Definition 6.1 We ^/z«// ca// f/ze elements ofS+ (u) the NESS of (M, 0, cj). 



Two basic properties of a NESS cj + are: 
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Proposition 6.2 (1) uj + e Si. 

(2) u+(a) > 0. 

Proof. The statements follow from the relations ( 16.591 ) and d3.7t . □ 

The assumptions (NESS1)-(NESS4) naturally apply to a family (M, <f>x, Wx, &x) except that in this case in 
(NESS3) one also requires that $^ € C(M). 

Remark 1. In the study of specific models it is an important and often very difficult problem to prove that 
is a singleton, namely that there exists u> + € S such that for all / € C(M), 

1 f* 

lim - / w s (/)ds = 

Remark 2. The regularity assumptions (NESS2) and (NESS3) are made for simplicity of presentations and can be 
relaxed, see Sections [6.4. 2l andl9lfor examples. 

Remark 3. The NESS property is related to the SRB property in dynamical systems, see HRu2l for more details. 

6.2 The Gallavotti- Cohen Fluctuation Theorem 

Let (M, <fi, to) be a lTRII system and let oj + E S+(w) be given. Let 

e t+ (a) =\ogu: + (e~ a ti°° ds ). 

Note that e t +{a) is a convex function of the parameter a and that et+(a) > —atu+(<j). 
We suppose: 

(GC) The Gallavotti-Cohen functional (GC-functional) 

e+(a) = lim \e t+ (a) 1 

t— >oo t 

exists for all a£l. 
If the GC-functional satisfies 

e+(a) = e+(l — a), 

for all a, we shall say that the GC-symmetry holds. e+(a) is a convex function and the GC-symmetry implies that 

e+(0) = e+(l) = 0, e+(a) < Ofora e [0,1], ande+(a) > Ofora ^ [0,1]. 

In comparison with the ES-symmetry, we remark that in general the relation et+(a) = et+(l — a) does not hold 
for finite t and that GC-symmetry may fail even in some very simple models (see Subsection 18.31 ). In contrast, 
whenever e(a) exists, the universally valid finite time ES-symmetry e t (a) = ej(l — a) implies that e(a) = 
e(l - a). 

Proposition 6.3 (1) Suppose that e+(a) is differentiable at a = 0. Then \(E3j\ holds and 

1 f* 

(<t)+ = Lu + (a) = lim - / u(a s )ds = -e,(0). 

t->oo t ,/ 
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If in addition the GC-symmetry holds, then {<r)+ = iff e+(a) = Ofor a G [0, 1]. For any regular sequence 



1 /■*» 

lim — / cr s (a;)ds = (cr)+, 



for uj + -a.e. x G M. If (o~)+ > 0, f/;en uj and uj + are mutually singular. 

(2) Suppose that e + (a) is differentiable for all a. Then the Large Deviation Principle holds for a w.r.t. (M, </>, oj + ) 
wif/i f/;e concave rate function I+{s) = mi a ^(as + e+(a)). T/'f/ie GC-symmetry holds, then 

I + (s) = s + I + (-s). 

The last relation is called the GC-symmetry for the rate function I + . 

Proof. The only part that requires a proof is the last statement in Part (1). Suppose that (cr)+ > and let 
uj + = V\ + V2, v\ <C uj, V2 -L oj, be the Radon-Nikodym decomposition of lu + w.r.t. uj. Since uj + G Si, 
Assumption |(C)| and the uniqueness of the Radon-Nikodym decomposition imply that v%, v-i G Si. If v\ is non- 
trivial, then Corollary 15 . 5 1 implies that (o)+ = 0, a contradiction. □ 

Consider a family (M, <j> x , Ux, X G M w , satisfying [(TT)][(T2)] and let lo x + G be given. Let 

9t+ (X,Y) = \oguj x+ (e-r-ti**'*'} . 

We suppose: 

(GGC) The Generalized Gallavotti-Cohen functional (GGC-functional) 

g+(X,Y)= lim Jst+(A:,r), 

exists for all X, Y, 
If the GGC-functional satisfies 

g+(X,Y)=g + (X,X-Y), 
for all X. y, we shall say that GGC-symmetry holds. 

Proposition 6.4 (1) Suppose that Y M> g+(X, Y) is differentiable atY — Q. Then 

1 t* 

($x)+ = lim - / w x (# Xs )ds, 

t^OO t J Q 

exists and = = — Vy.g+(X, F)|-k=o- Moreover, for any regular sequence t n , 

1 /•*» 

lim — / * Xs (2;)ds = 

for oj x +-a.e. x G M. 

(2) Suppose that Y i— > g + (X, Y") is differentiable for all Y . Then the Large Deviation Principle holds for Q x w.r.t. 
[M, (j) X , uj X +) with the concave rate function I x+ (s) = infy eR jv(Y • s + <7_|_(X, Y)). /ff/ie GGC-symmetry 
holds, then 

I x+ {s)=X-s + I x+ (-s). 
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(3) Suppose that g+(X, Y) is C 12 in a neighborhood of (0, 0) and that the GGC-symmetry holds. Then the kinetic 
transport coefficients are defined and satisfy the Onsager reciprocity relations Ljk = Lkj- 

(4) In addition to the assumptions of (3) suppose that for some e > 0, 

sup \\g t+ {^Y)\ <oo. (6.60) 
YeDe,t>i t 

Then the Fluctuation-Dissipation Theorem holds. 

The proof of this proposition is the same as the proof of Proposition l6.4l Apart from the last statement in (1) and 
(2), the conclusions of these two propositions are also identical. Note that gt+(0, Y) = gt(0, Y) and so the parts 
(4) of the two propositions are in fact identical (we included the statement for completeness). 

We shall say that a given ITRj system satisfies the Gallavotti-Cohen Fluctuation Theorem if the respective func- 
tional e_|_ (a)/g+ (X, Y) exist and are differentiable/C 1,2 and satisfy the GC/GGC-symmetry. It follows from 
Propositions 16.31 and 16.41 that the Gallavotti-Cohen Fluctuation Theorem is also an extension of the fluctuation- 
dissipation theorem to the far from equilibrium region. 

6.3 The resonance interpretation of e + (a) 

Let uj + be a NESS of (M, <f>, w) (in particular, we assume that (NESS1)-(NESS4) hold). For p € (-oo, oo], p ^ 
and / € C(M) let 

One easily shows that: 

Proposition 6.5 (1) uj + ([U p (t)f][U. p (t)g]) = u+(fg). 

(2) U p (t) extends to a strongly continuous group of bounded operators on L 2 (M,duj+) such that U*(t) — 
U-p(—t) and 

\\U p (t)\\ < e i*i" l -+, 

where m p+ = sup x6M | cr(ac) | / jjpj. Let L p+ be the generator of U p (t), U p (t) = e tLp+ . Then L* + = —L^ p+ , 
sp(L p +) C {z | \Re z\ < m p+ }, Dom(L p+ ) = Dom(i oc+ ), and for f e Domi p+ 

L p +f = L 00+ f H /. 

P 

We shall call the operator L p+ the NESS L p -Liouvillean. If a = — 1/p, then 

e t+ (a) = log(l,e tL "+l) + = log / e tL '+ldu+. 

JM 

Under suitable regularity condition this relation leads to the identification of e + (a) with a complex resonance of 
L p+ . With the obvious notational changes Propositions 15.161 and 15.171 apply to e + (a) and L p+ . 
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6.4 Examples 

6.4.1 The microcanonical ideal gas 

In the example of Section [3.6.1l it follows from ( 13.221 ). ( 13.23b that the unique NESS of the system is 

do, + = n^;-^)^- 

i=i 

We note that it is singular w.r.t. the reference measure u>. One immediately computes 

e +t (a) = -at(a} + , 

and observes that the GC-functional 

e+ (a) = -a(a) +) 

does not satisfy the GC-symmetry. Entropy production does not fluctuate in the NESS Accordingly, the rate 
function for its large deviation is 

f if s = (a)+, 
I — oo otherwise. 

The GC-fluctuation theorem fails in this model. 

6.4.2 The harmonic chain 

To compute the NESS in the example of Section l5.5.2l (in the thermodynamic limit) we note that u>xt is the centered 
Gaussian measure of covariance 

D t = (e- t£ *(/3/ l -fc(X))e- t£ )- 1 = h~^ a (fi - h-^e-^'^X^e-^h- 1 ^)- 1 ^ 1 ^. 

By Equ. ( 15.57b . this covariance converges strongly to the limit 

D+ = h- 1/2 W-(/3 - X^Wlh- 1 ' 2 . 

It follows that the system has a unique NESS uix+ which is Gaussian with covariance D + . In particular, one has 

gt+(X,Y) - -ilogdet (i + J D+e sC * 4>(Y)e sC As 
As for the functional gt(X, Y), one can compute the infinite time limit and get (see [JLTP]) 

g+(X,Y) = lim -a t+ {X,Y) = g(X,Y), 

t—>oo t 

i.e., the GGC-functional coincide with the GES-functional. In particular, it satisfies the GGC-symmetry and the 
Gallavotti-Cohen fluctuation theorem holds. Note that in this example Assumptions |(NESS2)| and |(NESS3)| do not 
hold. 
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7 The principle of regular entropic fluctuations 

The mathematical similarity between Propositions 15. 141 15. 151 on one side and Propositions 16.31 [6~4l on the other 
side is not accidental. The principal distinction is that the Evans-Searles symmetries are universal — they hold 
for any TRI dynamical systems for which the objects in question are defined. The mechanism behind Gallavotti- 
Cohen symmetries a priori could be model dependent and in general they may fail. Perhaps surprisingly, a careful 
look at all principal classes of models for which the symmetries have been rigorously established reveals that the 
respective functionals satisfy 

e(a)=e+(a), g(X,Y) = g+(X,Y). (7.61) 

This is true for the toy models (see Section©, for Hamiltonian open systems systems treated in BRT11 (see Section 
14.41 ) and for Anosov diffeomorphisms of compact manifolds BGCUlGC2l (see SectionfTTI). In each single case the 
strong ergodic properties of the model force the relations d7.6U . Note that the first relation in d7.61| ) holds iff the 
limits in the expression 

e + (a) = lim lim ^ log w u (e"%* <Tsds ), (7.62) 

t— >oo u— foo t 

can be interchanged and a similar remark applies to the second relation. This leads to a transparent mechanism 
for validity of GC-symmetries: The strong ergodicity (chaoticity) of the model forces the identities ( I7.6U and 
the universal ES-symmetries imply that the GC-symmetries hold. This leads to the principle of regular entropic 
fluctuations. 

Definition 7.1 We shall say that (M, <f>, uj, lu + ) satisfies the principle of regular entropic fluctuations if e(a) and 
e + (a) exist, are differentiable in a neighborhood of a — 0, and satisfy 

e(a) = e + (a), 

for all a. In the presence of control parameters, we say that family (M, (j)x, uJx, &x+, &x) satisfy the principle 
of regular fluctuations if g(X, Y) and g+(X, Y) exist, are C ' in a neighborhood of {Q, 0) and satisfy 

g{X,Y)= 9+ {X,Y). 

Combining Propositions 15. 141 15. 151 [6731 |6.4| one derives the implications of the proposed principle. The principle 
will naturally extend to quantum statistical mechanics. 

We emphasize that the principle of regular entropic fluctuations is ergodic in nature. It would be very interesting to 
exhibit examples of non-trivial physically relevant models for which the principle fails in the sense that e(a) and 
e+(a) exist and are differentiable, the GC-symmetry holds, but e(a) 7^ e + (a). 

If the principle of regular entropic fluctuations holds in a given model, then one may consider the Evans-Searles and 
Gallavotti-Cohen symmetries as mathematically equivalent. This does not mean that implications of Evans-Searles 
and Gallavotti-Cohen Fluctuation Theorems are identical. If (M, <j), uj) has non vanishing entropy production, then 
ui and uj + are mutually singular, and the corresponding Large Deviation Principles (LDPs) are very different 
statements. The principle only asserts that these LDPs hold with the same rate function and explains the origin of 
the Gallavotti-Cohen symmetry. 

Note that the Evans-Searles fluctuation theorem requires less regularity than the Gallavotti-Cohen fluctuation theo- 
rem and in particular does not require the existence of a NESS. This does not make the notion of NESS redundant — 
the fine studies of systems far from equilibrium are critically centered around the NESS. The situation is somewhat 
analogous to studies of phase transitions in spin systems — the pressure functional provides some information about 
the phase transitions but a much finer information is encoded in the set of equilibrium states at the critical temper- 
ature. 
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In Sections [10] and [TT] we shall examine the ergodic mechanism behind the principle of regular entropic fluctua- 
tions in the case of chaotic homeomorphisms of compact metric spaces and in particular Anosov diffeomorphisms. 
Starting with the works IGClllGCZ I these models have been a basic paradigm for any theory of entropic fluctua- 
tions. 



8 Toy models 

In this section we discuss a few additional examples for which the Evans-Searles and Gallavotti-Cohen functionals 
can be computed explicitly. 



8.1 Bernoulli shift 

Let M = {0, 1} Z . The elements of M are sequences x — (xj)j € %, Xj S {0, 1}. Let v p be the usual Bernoulli 
measure on {0, 1}, v p ({l}) = p, ^({0}) — I — p. Let 

w = ( (8) v p J ( ® v * J 

\n— — oo / \n— 1 / 

where p, q S]0, 1 [ are given. The dynamics </> is the left shift 

<f>{x)j = Xj+i. 

The Radon-Nikodym derivative of tJi = uj o w.r.t. ui is 

A , (x) = Vpi{xo}) 

and the entropy production observable is 

Au f lo §? if*o = l, 

a(x) =logA Ul \ u (x) =log— A (x ) = < 

d »P llosfEf ifx = 0. 



Note that for all n > 1, 



w(er n ) = qlog- + (1 - q)log] = -Eni,(v g \i/ p ) > 0, 

P 1 ~P 

so that the system is entropy producing, 

(a)+ = -Ent(v q \v p ) > 0, 

iff p q. Similarly, for all n > 1, 

ilogw(e a £"=i'-i) =log[p 1 - a g a + (l-p) 1 - a (l-g) a ] = Ent Q (^ p ), 
and so the ES -functional exists and is given by 

e(a) = log [p 1 - a q a + (1 - p) 1 ^! - q) Q ] . 
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Assuming p ^ q, the ES-symmetry e(a) = e(l — a) holds iff g = 1 — p. In particular, if q ^ 1 — p, the model is 
not time reversal invariant. If q = 1 — p, then the time reversal is r d(x)j = 1 — x-j. 

Let 

M is a compact metric space and for any / E C(M), 

lim uj n (f) = 

n—t-oo 

Hence, cj+ is the NESS of (M, 0, cj). Moreover, for any n > 1, 

i log W+ (e- a ^"=o = log [tf 1 " V + (1 - g) 1 "^ - p)«] , 
and so the GC-functional is 

e+(a) = e(l — a). 

The GC-symmetry holds iff the ES-symmetry does, and in this case the GC-functional coincide with the ES- 
functional. 

If p = 1 — q, then one can consider X = p — 1/2 as a control parameter. The respective Fluctuation-Dissipation 
Theorem follows easily IShll . 



8.2 Baker transformation 

Let M = [0, 1 [x [0, 1 [, duo = dxdy and 



(x/p,yq) ifxe [0,p[, 

((x-p)/q,py + q) ifx£[p, 1[, 



where p, q e]0, 1[. This model has been studied in iDGTl . see also (Do] |TGl EB . Define a map T : {0, 1} Z -> M 

by T{a 3 -}j e z = 0, v) where 



S = P ( a + J2 OnP 71 -^' 0< <? E? =° °* ) , y = 9 U-i + J2 a-n-iq 71 -^ a — 'q^ a — 1 

n=l / \ n=l 



The map T is a mod isomorphism between the Bernoulli shift considered in the previous section and the Baker 
transformation. Hence, the ES-functional e(a) of these two models are the same. The non-equilibrium steady 
states of the two models are also related by the map T and their GC-functionals are equal |Sh|. 



8.3 Dilation on a half-line 



Let M = [0, oo] (the compactified positive half-line), (f>t(x) — e yt x where 7 6 M, and 

2 dx 

duj = — 5- . 

n 1 + x z 

The map d(x) = x^ 1 is a time-reversal of (M, <j>j,uS). The reference state oj is invariant under <fio and for 7 7^ 
the system has a unique NESS 



LJ 1+ — 




if 7 < 0, 
if 7 > 0. 
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One has 

A^|^(x) = e~ 7 * _^ cr 7 (a;) = -7-— — ^ 

i 1 + e 1 + ar 

from which it follows that (c 7 ) + = cj 7+ (er 7 ) = |7|. One easily computes the ES-functional 

e(a) = 



-a\-y\ if a < 1/2, 

-(l-a)|7| ifa>l/2, 



while the GC-functional is e+(a) = — a|'y|. Hence, e(a) ^ e+(a) for a > 1/2, the ES-symmetry holds, but GC 
does not. 

Consider 7 as a control parameter. The associated flux observable $ = (x 2 — l)/(a; 2 + 1) does not depend on 7, 
u;(<£>) = and w 7+ (<l>) = sigrry. The kinetic transport coefficient 



d-f w 7 +(*)| 



7=0 ' 



is not well defined. On the other hand, w($ t ) = th(7i/2) is a real analytic function of 7, the finite time kinetic 
transport coefficient is well defined, the finite time Green-Kubo formula holds, and 



d 7 



Note that, as in the example of Section 15.5.11 L t diverges as t — > 00. The reader may have noticed the similarity 
between this toy example and the microcanonical ideal gas of Section [3.6.1l (see also IICGII for a related example). 

By slightly modifying this example we can illustrate another point. Consider (M, </>_ 7 2 , uj) and let again 7 be the 
control parameter. Then 

= So for all 7/O, and the functions (c 7 )+ = -f 2 , ($ 7 )+ = 7 are entire analytic. The kinetic transport 
coefficient is equal to 1 and, since $0 = 0, the Green-Kubo formula fails. In this example the finite time linear 
response theory holds with L t — 0, 7 1— > (3> 7 )t is real analytic, 7 1— > ( < 1 ) 7 )+ is entire analytic, but the limit and 
derivative in the expression 9 7 (lim t _ s . 00 ($ 7 ) t )| 7= o cannot be interchanged. 



9 Gaussian dynamical systems 

This class of dynamical systems is treated in detail in the forthcoming article IIJLTPI and for reason of space we 
shall be brief (in particular we will omit all the proofs). The thermally driven harmonic chain of Section IB. 5. 21 is 
an example of a Gaussian dynamical system. 

Let r be a countably infinite set and 

M = R r = {x = {x n ) ner I x n G K}. 
A sequence I — {l n } n< zr of strictly positive real numbers such that X) n er — 1 defines a metric 

|*^n Vn I 

ner 



= 5 1 + 



on M. Equipped with <i, ill is a complete separable metric space. Its Borel cr-algebra J- is generated by the 
cylinders 

C(B- m, . . . , n k ) = {x G M I (x ni , . . . , x n J G B}, (9.63) 
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for k > 0, m, . . . , n/- e T and Borel sets B cR k . 

We denote by £|(T) C M (respectively M; C M) the real Hilbert space with the inner product (x,y) — 
Sner x nUn (respectively (x, y)i = J2 n er InXnVn)- ^«(r) is dense in A// and Mi is dense in M. All the measures 
on (M, F) we will consider in this example will be supported on Mi. We denote by A nm = (S n , AS m ) the matrix 
elements of a linear operator A on £^(T) w.r.t. its standard basis {S n } ne r- 

Let £ be a bounded linear operator on £j|(r) which has a continuous extension to M;. For x € M and i € K we 
set 

!e t£ a; if x £ M h 
(9.64) 
a; if x ^ Mi. 

0* is a group of automorphisms of (M, J 7 ) describing the time evolution. Note that the map (t, a;) n- is 
measurable and so Assumptions |(Fl)|(F2)| of Section l2~3l hold for (M, (j>). 

Let 13 be a bounded, strictly positive operator on £^(T). The centered Gaussian measure on (M, J 7 ) of covariance 
D is the unique measure cj specified by its value on cylinders 



which shows that uj(Mi) = 1, 

Our starting point is the dynamical system (M, </>, w). Lu t — lo o <\r l is a Gaussian measure of covariance 



£> t is a bounded strictly positive operator on £j|(r) and U)t(Mi) — 1 for all t. Denote by T the real vector space 
of all trace class operators on The trace norm ||T||i = tr((T*T) 1 / 2 ) turns T into a Banach space. By the 

Feldman-Hajek-Shale theorem uj t and u) are equivalent iff D^ 1 — D^ 1 E T. We shall assume more: 

(Gl) The map K 3 ^ D^ 1 - L^ 1 e T is differentiable. 

It follows that 




where D c — [D ninj ]i<i : j<k. For any finite subset Acf one has 




D t = e tL De' 




^- 1 )|t=o 



is trace class. Set 



a(x) = (x,<;x) 



tr(Z^)- 



(9.65) 



Proposition 9.1 Suppose that (Gl) holds. Then: 



(1) cr G L 1 (Af,dw) anc/i i— > er t is strongly continuous in ^(Mjdiv). 

(2) ^ t |cj = J cr- s ds ami i h-> e^^ti" ;s strongly C 1 in L 1 (M, da;). 

(3) w t (cr) = tr(<r(D( — D)) and in particular u{a) = 0. 



(4) Ent(w t |w) - - /J tr(?(D 8 - £>))ds. 
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This proposition implies that Assumption |(E 1 )| holds and that a is the entropy production observable of (M, <fi, to). 
In some examples only finitely many matrix elements q nm are non-zero and in this case a is continuous, a is 
bounded only in the trivial case q — and so Assumption |(E2)| is not satisfied. 

Our next assumptions are: 

(G2) For some constants rn± and all f € M, < m_ < D t < m + < oo. 
(G3) The strong limit 

s — lim D t = D + 

exists. 

Clearly, m_ < D + < m + . Let oj + be the Gaussian measure on (M, J-) with covariance D + . 
Proposition 9.2 Suppose that (G1)-(G3) hold. Then: 

(1) For all f e C R (M), 

lim u) t (f) = u+(f). 

(2) a E L 1 (M,duj+) and 

(a)+ = lim w t (cr) = tr(?(D+ - D)) = w+(tr). 

We shall call w + the NESS of (M, 0, w). Note that assumptions (NESS2) and (NESS3) of SectionOdo not hold 

for (M,<f>,w). 

Finally, we assume the existence of time reversal in the following form: 

(G4) The exists a linear involution d on ^g(F) such that fiC = -Cd and dD = D-d. 

This assumption implies that D- t = flDtd also satisfies (G2) and (G3) and 

L>_ = s - lim D t = ■dD+'d. 

t— > — oo 

Moreover, $q — — q$ and ti(Dq) = 0. 
Since a is unbounded, 

e t (a)=logw(e 0, /<f ,r - d *) ) 
is a priori finite only for a E [0, 1]. Note that et(a) is real analytic on ]0, 1[. (G4) implies that 

et(a) =togw(e- a /o 

and that e t (a) — e t (l ~ a). Set S — m-/{m + — m_). 

Proposition 9.3 Suppose that (G1)-(G4) hold. Then: 

1. et{a) is finite and real analytic on the interval ] — 8, 1 + S[. 

2. D a = ((1 — a)D^ + aDZ 1 ) is a real analytic, bounded operator valued function of a on this interval. 
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With these preliminaries, the Evans-Searles Fluctuation Theorem holds in the following form. 
Theorem 9.4 Suppose that (G1)-(G4) hold. Then: 
(1) Forae] - 5,1 + 5[, 



1 f a 
t e ^ = -J 



e(a) = lim -e t (a) = - / tr(<;D 7 )dj. (9.66) 



t— >oo 

The ES-functional e(a) is real analytic and convex on the interval ] — 1 + S[, satisfies the ES-symmetry and 
e?(0) = -(*)+. 

(2) Ift n is a regular sequence, then 

1 /■*» 

hm — / <j s (x)ds = (cr)+, 

n-S-oo t n J 

foruj-a.e. x. 

(3) The Large Deviation Principle holds in the following form. The function 

7 ( S ) = ^ (as + e(a)). 

is a concave with values in [—00, 0], I(s) = iff s = (cr) + and I(s) = s + I(—s). Moreover, there is e > 
such that for any open interval J c] — (<x)+ — e, (cr)+ + e[, 



lim - log co [ i x 



- f a s (x)ds G jl ] =sup/(s). 



We now turn to the Gallavotti-Cohen Fluctuation Theorem. Let 

e t+ (a) =log W+ (e- a/ ° CT8ds ). 
A priori e t+ (a) might not be finite for any a. 

Proposition 9.5 Suppose that (G1)-(G4) hold. Then: 

(1) e t+ (a) is real analytic on the interval ] — 5, 1 + 5[ and for any a in this interval, 

1 f a 

e+(a) = lim -e t+ (a) = — / tr(^D 7 )d7 = e(a). (9.67) 

t->oc t ,/ 
/« particular, (M, (j),uj,Lj + ) has regular entropic fluctuations. 

(2) /fi„ is a regular sequence, then 

lim -!- / <7 s (x)ds = (cr)+, 

for ui+-a.e. x. 

(3) 77ze Large Deviation Principle holds for a and (M, <j>,w+), i.e., for some e > one/ for any open interval 

J C] - - e, + e[, 



Um llog W+ ( >r 



7 / cr s (a;)ds e jl ) =sup/(s). 
1 Jo J / sgj 
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The proofs of all the results described in this section can be found in [JLTP| and here we will only sketch the 
computations leading to the Formulas ( 19.661 1 and j9.67t , 



Using 



d a e t {a) 



( J *cr s ds e" a So °>> d A 



uj( e - a So^ ds ) 

and the fact that e~ a So CTs dw /uj(e~ a So is a Gaussian measure of covariance [(1 — ajD^ 1 + aDZ]\ , we 
get 



e t (a) 



= -t 



J J tr(e fl£ V £ [(l- 7 )^- 1 +7^It] _1 )d S d 7 
£ J tr [(1 - ^D- 1 + _1 ) dsd 7 

(l-7) J Dr s 1 +7^: f 1 (1 _ s) l jd S d 7 , 



o -'o 



tr c 



and so 



lim — et(a) = — lim 
Regarding ( 19.671 ). we have 



a f l 



JO 



tr ? 



S a e t+ (a) 



(1 - 7 )D^ + 7^Z 4 1 (1 _ J d,sd 7 = -J tr( ?J D 7 )d 7 . 
w+ ([/o<r.d«] e- a /o*- d ») 



a; + (e- a ^o^ ds ) 

where e _Q ^o CTs dw + /u;(e~% <T a d;i ) is a Gaussian measure with covariance [-D+ 1 + <xDZt — a!? -1 ] . Pro- 



ceeding as before, we get 



e t+ (a) = -t 



tr <r 



^; 1 + 7^:! (1 _ S) -7A" S 1 



dsd7, 



and hence, 



lim — e t +(a) = — lim 

t— T-OO £ t— >00 



JO 



tr ? 



^; 1 +7^: t 1 (1 _ s) -7^ 1 



dsd7 



tr(cD 7 )d7. 



We finish with several remarks. 



Regarding the resonance interpretation of e(a) and e+(a), since a is unbounded the study of Liouvilleans and 
their resolvents requires some care. Regarding generalized functionals and symmetries, if Cx and Dx depend on 
control parameters X, then under mild additional regularity assumptions one can compute g(X, Y), g + (X, Y) and 
prove the Fluctuation Dissipation Theorem. Again, g(X, Y) = g+(X, Y) and the principle of regular fluctuations 
holds. 



10 Homeomorphisms of compact metric spaces 

Let [M, d) be a compact metric space and <j> : M — > M a homeomorphism. In this section we consider the 
discrete time dynamical system on M generated by <f>. We shall show, using the thermodynamic formalism, that 
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if (j) is sufficiently "chaotic", then the ES and GC Fluctuation Theorems, the Fluctuation-Dissipation theorem and 
the principle of regular entropic fluctuations hold for a large class of reference states ui. Our treatment generalizes 
HMVI [Mali. Some of the results presented here extend also to flows (see Remark 4 on Anosov flows at the end 
Section [TQ. 

The material in this section is organized as follows. Subsection 110.11 is a brief review of some basic aspects 
of topological dynamics (see [Wal | for a more detailed introduction). Subsection 110.21 deals more specifically 
with two classes of "chaotic" topological dynamics: expansive homeomorphisms with specification (see IKHIO 
and Smale spaces (see llRull ). The reader familiar with these topics can skip Subsections 1 1 0. 1 1 1 0.21 and proceed 
directly to Subsection 110. 31 where, adopting the point of view of IMVl lMaTl . we discuss entropy production and 
fluctuation theorems for chaotic homeomorphisms. We show in Subsection 110.41 how these results relate with 
the general approach, based on the notion of reference state, advocated in this work. Finally, we discuss, as an 
example, the simple case of a topological Markov chain in Subsection |10.5l 

10.1 Topological dynamics 
10.1.1 Entropy and pressure 

Let v E Si and let £ = (C\,C2, ■ ■ ■ ,C r ), Cj E J 7 , be a finite measurable partition of M, By a standard 
subadditivity argument, the limit 

MM) = - lim I ]T u(r\c j Jn---ncf>-\c j J)io g u(4>- 1 (C h )n---nr k (C jh )), 

j'ii...,i*6{l, —,r} 

exists. The Kolmogorov-Sinai entropy of <j> w.r.t. v G Sj is given by 

h v {<t>) = supfc„(&0- (10-68) 

i 

The map v M> h v ((f)) is affine, that is, 

M+(1-A) M (0) = *h v {<t>) + (1 - 

for any v, \i E Si and A <G [0, 1]. 
For x E M and e > we denote by 

B± n (x, e) = {yE M \ max d{4> ±j {x), cj> ±j (y)) < e}, 

0<j<n— 1 

the Bowen ball of order ±n . The following result is known as the Brin-Katok local entropy formula |BK| and is 
a topological version of the Shannon-McMillan-Breiman theorem. 

Theorem 10.1 Suppose that v E Sj is non-atomic and that h v (<j)) is finite. Then for v-a.e. x E M, 
limlimsup log v(B n (x, e)) = lim lim inf \ogv(B n (x, e)) = h v {x). 

t-s-0 r n.oo n e-s-0 n->oo n 

The function h u (x) is (^-invariant and v(h v ) = h v (<p). If v is ergodic, then h v (x) — h v (<f>) for v-a.e. x. 
For any function <p E C(M) we set 

n-l 

s n ip(x) = y^(po^(x). 

3=0 
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A set E C M is said to be (n, e)-separated if for every x, y £ E, x 7^ y, we have &/ ^ B n (x, e). For G Cr(M) 
let _ 

Z n (^e) = sup^e s "^, 



E xeE 

where the supremum is taken over all (n, e)-separated subsets of M, The limit 

P(<p) = lim limsup- log Z n (<p,e), (10.69) 

e— ^0 n— ¥00 Tl 

exists and P(<p) is called the topological pressure of with respect to the potential (p. Alternative representations 
of the pressure are 

P{ip) = lim limsup -log V e s " v{x) = limliminf -log V e s " v( - x) , (10.70) 

where the E n>e are arbitrary maximal (n, e)-separated sets. The pressure also satisfies the variational principle 

P(<p) = sup (v(<p) + K(4>)) . (10.71) 

The special case P(0) = sup^^^ /ij/(</>) is the topological entropy of </>. An immediate consequence of the 
variational principle is that P(p) depends only on the topology of M and not on the choice of metric d. Another 
consequence is that the map Cr(M) 3 ip -P(</?) is convex and that either P((p) = +00 for all ip, or P(ip) is 
finite for all ip. In what follows we assume that P{p) is finite for all p. Occasionally we shall use an additive 
normalization (p = p — P(<p) which ensures that P(p) = 0. 

10.1.2 Potentials and equilibrium states 

An invariant measure v is called an equilibrium state for the potential p if the supremum in ( 110.71b is realized at 
v. In the case ip = the equilibrium states are called measures of maximal entropy. In general, equilibrium states 
do not necessarily exist, and if they exist they are not necessarily unique. The set S cq (ip) of all equilibrium states 
for p is obviously convex. If S cq (p) is singleton, we denote by v v the unique equilibrium state for ip. 

Theorem 10.2 Suppose that the entropy map 

Si3v^h v {4>), (10.72) 

is upper-semicontinuous. Then: 

(1) For all p £ Cr(M) the set S cq (p) is non-empty and compact. A measure v is an extreme point of S cq (p) iff 
v is <j)-ergodic. 



(2) For a dense set of p in Cr(M) the set S cq (p) is a singleton 

cq 



(3) For p £ Cr(M), the map R N 3 Y i-> P(tp + Y ■ f) is differentiable at Ofor all f £ Cr(M) n ijfS cq {p) is a 



singleton and in this case 



W Y P(p + Y-f) 



Y=0 
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The proofs of (1), (2) and (3) for N = 1 can be found in [Wal |. For a proof of (3) for N > 1 see | Je|. 

For a given ip £ Cr(M) we denote by S v the collection of all v € <S such that for all x € M, n > and 
sufficiently small e > 0, 

Cn(e)- 1 < v{B n {x,e))e- s ^ < C„(e), (10.73) 

where C„(e) > satisfies 

lim -logC*„(e) = 0. (10.74) 

n— too Ji 

The main property of the class S v is the following basic result of Kifer lKH which relates it to the large deviation 
formalism trough the topological pressure. 

Proposition 10.3 For v £ S v and ip € Cr(M), 

lim - log v(e s ^) = P(<p + tp). (10.75) 

n— foo fi 

Proof. Let E n ^ t be arbitrary maximal (n, e)-separated sets. Note that for two different x, y € E n _ e one has 

B n (x, e/2) n B n {y, e/2) = and by maximality [J xeEn ^B n {x, e) = M. It follows that 

where 1 s denotes the indicator function of the set B. Setting 5 t — sup x yeM d / x y \ <e \ip(x) — ip(y)\, we obtain 
E v{B n {x,e/2))c s ^^- nS ' < v(e s ^) < £ v{B n {x,e))e s ^^ +n5 < . 

Combining these estimates with ( 110.731 ) we get 

Cnie/2)- 1 £ e 5 "^)^-" 5 ' < ^(e 5 "^) < C„(e) £ e 5 "^+^^+"S 

and so 

lim sup - log E e s "^ + ^ (a:) > lim sup - ]ogv(e 3 ^) - <5 e , 
liminf - log V e^ 4 ^^ < liminf i logz/(e s "^) + <5 e . 

n— ^oo tt, ^ — ' n— >oo 77, 

Since <5 e | as e | 0, the statement now follows from Equ. ( 110.701 ). □ 

Corollary 10.4 Suppose that the entropy map M0.72\ is upper-semicontinuous and that S cq (cp) is a singleton. 
Then, for all v £ S^, 

lim -S n f(x) = v v (f), 

n— too fi 

for all f €E Cjr(.M) and v-a.e. x G Af. /n particular, is the unique NESS of the system (M , 0, 
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Proof. By Proposition [TU3J the generating function 

lim -logv(e- aS " f ) = PUp- af), 

exists for all a E M.. By Theorem l 10.21 it is differentiable at a = so the statement follows from the Gartner-Ellis 
Theorem (Proposition |5.1 11 (1)). □ 

As an immediate consequence of the Katok-Brin local entropy formula (Theorem llO.il ) we have: 
Proposition 10.5 Suppose that v € Si n S v is non-atomic. Then for small enough e and v-a.e. x € M, 

h v (x) = — lim —logf(B n (x,e)) — — lim —S n <p(x), 

n— too fi n— too TL 

and in particular f((p) + h v (<p) = 0. 

Corollary 10.6 If S v is non-empty, then P(ip) = 0. Moreover, if v E Si C\ S v is non-atomic, then v E S eq (ip). 

Proof. Setting if> = in ( 110.75b we get P{<p) = 0. This fact and Proposition ! 10.5 l imply the second statement. □ 

We shall say that the potentials tp and %p are physically equivalent, denoted ip ~ tp, if 

lim — sup liSn^fac) — S n tp(x)\ = 0. 
n->oo n xeM 

This clearly defines an equivalence relation on Cr(M). The following facts are easy to prove: 

Proposition 10.7 (1) If tpi ~ ant/ 952 ~ "02, a<pi + b(f2 ~ a^i + b^ifor all a, 6 e K. 

(2) (p ip implies P(<p) = P(tp) and u(cp) = u(tp)for all v E 5/. /n particular, S cq (ip) — S e<l (ip), 

(3) S^^S-^iffip^ ip. 

(4) Either S v C\ S^ = $ or S v = S^. 

(5) A^ofe f/za? <S V depends on the choice of metric d. If d is a metric equivalent to d (i.e., C~ 1 d < d < Cdj, then 
S v ,d = S^j. 

10.2 Chaotic homeomorphisms 
10.2.1 Expansiveness and specification 

The non-triviality of S v can be deduced from suitable "chaoticity" assumptions on <fi. A homeomorphism <p is 
called expansive if 

(ESI) There exists r > such that if d((p n (x), (p n (y)) < r for all n E Z then a; = y. 

r is called expansive constant of <p. If <p> is expansive then h v (<f>) = h v {<p, £) for any measurable partition £ such that 
diam(£) = sup{(i(x, y) \ x, y G C, C E £} < r and the entropy map 5/ 3 v n- h u (<p) is upper-semicontinuous 
(see Proposition 6.5 in [Rul |). 

cj) is called an expansive homeomorphism with specification if in addition to (ESI) 
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(ES2) For each S > 0, there exists an integer p(5) > such that the following holds: if a < b are 
integers, , . ,I n finite intervals of Z contained in {a, . . . , b} with dist(Ij,Ik) > p(S) for j ^ k, 
and xi,...,x n <E M, then there is x g M such that b ~ a +p( 5 ) (x) = x, and 

for k € Ii, i = 1, . . . , 7i. 



A potential i/3 is called regular if for all sufficiently small e > there exists C e > such that for all x £ M, all 
n > and all y G £?„ (x. e), 

\S n (p(x) - S n (p(y)\ < C t . 

If ip is regular so is = ip — P((p). 

Expansive homeomorphisms with specification were introduced by Bowen in [Bol ] and have been much studied 
since then. We recall the following classical result of Bowen HBolU (see also |KH|). 

Theorem 10.8 Suppose that 4> is an expansive homeomorphism with specification and ip a regular potential. Then 
S cq (tp) is a singleton and v v € S^. 



An expansive homeomorphism with specification cf> has a rich set of periodic points which completely determine 
the equilibrium state v v of a regular potential ip. Indeed, 

Mf) = lim ~V~T~\ E e^>/(.T), (10.76) 
xeFix(*») 

for all / € C(M), where 

z n (<p) = eS,Mx) > 

a;GFix(0™) 

and Fix(0") = {x € M \ (j> n (x) = x}, the set of periodic points of <fr of period n (see IIBoll ). Moreover, the 
pressure of <p is given by 

PQp) = lim -logZ n (<p), (10.77) 

n— foo ti 

(see e.g. [KHj, Proposition 20.3.3). These two approximation results lead to the following characterization of 
physical equivalence. 



Proposition 10.9 Suppose that <fi is an expansive homeomorphism with specification and that (p, -0 are regular 
potentials. Then the following statements are equivalent: 



(1) <p ~ 0. 

(2) For all n and all x G Fix(</> n ), S n ip(x) = S n ip{x). 



Proof. For x £ Fix(0") and k € N, one has Skn{<p — = kS n (tp — ilf){x) and hence 

\S n (<p - ip)(x)\ <n^— sup \Skn(<P ~ 
kn xeM 

Letting k — >• oo shows that (1) (2). Suppose (2) holds, then Equ. ( 110.77b implies P((p) = P(ip) and Equ. 
(110.76b v v = v^. Theorem 1 1 . 8 1 now implies that g n <S>^ and Part (5) of Proposition ! 10. 7l vields 5^ = S^. 

By Part (4) of the same proposition we have <p ~ -0 and, since -P(<ys) = P{i/j), we conclude that ~ 0. □ 
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10.2.2 Smale spaces 

One can say more under a stronger "chaoticity" assumption. (M, cf>) is called Smale space if the following holds: 

(51) For some e > 0, there exists a continuous map 

[•,•]: {{x,y) £Mx M\d(x,y) < e} -» M, 

such that [x, x] = x, [[x, y],z] — [x, z), [x, [y, z]] — [x, z] and 4>([x, y]) = [4>{x), (p(y)] whenever both 
sides of these identities are defined. 

(52) For some <5 > 0, 0<A<1 and all n £ N one has 

oo 

d(tf> n (y), 4> n (z)) < \ n d(y, z) if y,z G f| B m (x, S), 

m=l 
oo 

d{4>- n {y),4>- n {z)) < \ n d(y, z) if y t z e f| B_ m (x, 5). 

m—l 

A Smale space is called regular if there exists C > such that d(x, [x,y]) < Cd(x,y). A Smale space is 
topologically + transitive if there exists x € Al such that the set {<fi n (x)\n > 0} is dense in M. (M, </>) is a 
topologically + transitive iff for any open sets U, V and any N > there exists n > N such that <fi n (U) flV ^(5. 
A Smale space is topologically mixing if, for any open sets [/, V there exists iV > such that, for all n > N, 
4> n (U)nV 0. If (M, <fi) is a topologically + transitive Smale space then <fi is an expansive homeomorphism with 
specification. Note that topologically mixing =>■ topologically + transitive. 

Smale spaces can be studied using powerful tools of symbolic dynamics and are very well understood. We will 
recall some classical result (see Chapter 7 in MRul l and in particular Corollaries 7.10 and 7.12). For a g]0, 1[ we 
denote by C^ C (M) the real/complex vector space of all Holder continuous functions with exponent a, i.e., all 
/ e C R/C (M) such that, for some C > and all x, y € M, \f(x) - f(y)\ < C d(x, y) a . The norm 

||/||a = 8up l/( ; ) ~ff )l + supl/( 3 ;)l, 
X1 ty d(x,y) a x 

turns Cg/ C (M) into a real/complex Banach space. 

Theorem 10.10 Let (M, </>) be a topologically + transitive Smale space and let a g]0, 1[ be given. 

(1) The map C^(M) 3 ip t— > P(ip) is real analytic. 

(2) For any ip £ C^(M), S eq (tp) is singleton and v v £ S^. 

(3) Ifip, ip £ (M), then v v = iffip andip are homologous, i.e., ip = ip+c+fo<fi—f, where c = P(ip)—P(ip) 
and f £ Ck(M) is unique up to an additive constant. If(M, (f) is regular, then f £ C'^(M). 

(4) Ifip,ip £ C£(M), then ip ~ ip iff ip and ip are homologous and P(<p) = P(ip). 
Suppose that (M, <j)) is topologically mixing. Then 

(5) Iff, f,g£ C^t(M), then for some A, B > and all n £ 7L, 

IM5/n)-Ms)M/)l<Ae- B w. 
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(6) Suppose that <p, /W, ■ • • , /W € C|(M). Then the Central Limit Theorem holds for f = (/W, ■■■ , /W) 
w.r.t. (M, (j), v v ) with covariance matrix 

Moreover, if (M, 4>) is regular, then D^k > unless is homologous to 0. 

(7) Letip,?p € C§t(M). Consider the transfer operator 

on Cr(M) and let 

oo 

K(z) = yV"%(f7£l), Rez > Bup|^(x)|. 

n— 

Then for some e > the function IZ(z) has a meromorphic continuation to the half-plane 

Rez> P{tp + ip) -P{f) -e, 
and its only singularity is a simple pole at P(tp + ip) — P((p). 

(8) Let ip € (M). Then there exists e > and C e > rac/i that, for all ip £ C£(M) with \\ip\\a < e. 

sup -| log ^(e s "^) | <C e . 

For (l)-(6) see |Rul |. (7) and (8) are implicit in BRulllBafl and are easily established using the well-known spectral 
properties of Ruelle transfer operators. Note that (8) and Proposition I5T91 yield CLT for , . . . , / (Ar) G C^(M) 
w.r.t. (M, 0,1/^). 

10.3 Entropy production 

In this section we suppose that is lTRII with a continuous time reversal fl. Note that the map v ^ v o d preserves 
Si. One easily checks, using the definition ( 110.68b . that 

/w(0) = M<r X ) = hu(<f>). (10.78) 

It follows that, for any ip £ Cr(M), 

P(^s) = sup {v(ip) + h v {4>)) = sup o #(<^) + h vo $((j>)) 

= sup (1/(^3 o ft) + h v {<t>)) = PQp o i?). (10.79) 

v&Sj 

To each potential <p £ Cr(M) we associate the function 

= <fi - <f O <&. (10.80) 

As we shall see in the next section, a v is closely related to the entropy production observable a of the dynamical 
system (M, tj>, w) (as defined in Equ. (13 . 1 5b ) for uj E S v . In this section, we investigate the intrinsic properties of 
cL and its fluctuations. 
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Proposition 10.11 (1) a v o i9 = —<r v (compare with Proposition ^. 7\ . 

(2) Iftp ~ ^ 5 V = S^f,), then a v ~ Sty. 

(3) Suppose that S eq (ip) is a singleton. Then = iff o d — v v . 
In the remaining statements we assume that c £ 5/ fl 5^ is non-atomic. 

(4) For e small enough 



1 K^n(a:,6)) 1 
lim — log ; — t— — — — = lim — 

n->oo n vyB n (v o <p n L (x),e)) n-n»n 

/or v-a.e. x £ ill. 7/V is ergodic, then 



lim - log — = ,(5,), (10.81) 

Hoc n v(B n (v O (f> n L (x),€)) 

for v-a.e. x £ M. 

(5) v(jj v ) > 0. 

(6) Suppose that v is ergodic and that v and v o ■& are equivalent measures. Then f(<7„) = 0. 

Remark 1. Apart from (1), all statements of the previous proposition hold with the same proofs if one replaces -d 
by d o (f) k with an arbitrary k £ Z. Our choice of k — differs from the one in IIMVL A different choice would 
not affect any result in this and the next section. 

Remark 2. By part (4), the observable <7„ associated to a potential ip £ Cr(M) quantifies the Brin-Katok local 
entropy produced by changing the reference point from x to ■& o <j) n ^ 1 (x), i.e., by reversing the orbit of x. 

Proof. (l)-(2) are obvious. 

(3) v v (cr v ) = is equivalent to Vtp(ip) — is v (ip o $) which, by the variational principle and Equ. ( 110.781 ). is 
equivalent to 

P(tp) = v v (ip) + h Vv (4>) = v v o i%) + Knotty). 
Hence o ■& is also an equilibrium state for ip and the uniqueness implies v v = v v o -d. 

(4) follows easily from the conditions dl0.73l )-( fl0.74| i and Birkhoff ergodic theorem. If v € Si P\S V is non-atomic, 
then by Corollary |10.6l v £ iS eq ((/?). Since v o ■& £ Si, the variational principle and Equ. (110.781 1 lead to 

v{tp) + K{4>) = P{tp) >vo 0(<p) + h yoi) {4>) = v{Lp o <&) + h y {4>), 

and (5) follows. 

(6) The Brin-Katok formula and (110.811 imply that v(a v ) = h u (<p) - K(4>) = 0. □ 
The fluctuations of the observable a v in the states v £ S v are described in our next result. 
Proposition 10.12 (1) For all v £ S v the functional 

l9a4 e v (ot) — lim - log^(e- aS "^), 

n— j-oo 71 

exists and is given by e ip (a) = P(cp — aa^). 
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(2) The symmetry e v (a) = e v (l — a) holds. 

(3) Ifip ~ ip, then e„(a) = e^(a). 

(4) Suppose that the entropy map Sj 3 v t— > h v (<fi) is upper-semicontinuous and that S cq (ip — aa v ) is singleton 
for all a € M ana? Zef G S v . Then the Large Deviation Principle holds for a v w.r.t. [M, (fi, v) with the 
concave rate function I v (s) = sup QgR (as + e v (a)) which satisfies the relation I v (s) = s + I v (—s). 

Proof. (1) follows from Proposition ! 10. 31 Writing ip — aa v = (1 — a)ip + a(p o (2) follows immediately from 
Equ. dl0.79fc . (3) is a direct consequence of (1) and Proposition ! 10.71 (4) follows from Theorem ! 10.21 Proposition 
110.31 and the Gartner-Ellis Theorem (Proposition !5.1 U .D 

One can introduce control parameters in the above framework and discuss the generalized symmetry and linear 
response theory. Consider a map R N 3 X H> (fix, where each (fix is a homeomorphism of M with continuous 
time reversal dx, and a map R N 3 X H> ipx € Cr(M). Let ax — 5>x- We shall assume that there exists 
<& x = {<&x\ ■■■ , with g Cr(M), such that 

ax~X-<*>x, (10.82) 

and <&x ° $X ~ &x- Finally, denote by Six the set of </>x -invariant states and by Px the pressure functional 
for the map (fix ■ 

Proposition 10.13 (1) For all v € S Vx , the functional 

R N 3Y ^ q{X,Y) = lim -\ogv(e- s " {Y -* x) ), 

exists and is given by q{X, Y) = Px(fx — Y ■ &x)- 
(2) The symmetry g(X, Y) = g(X, X — Y) holds. 

In the remaining statements we assume that the entropy maps Six 3 v h u (<fix) ore upper-semicontinuous. 



(3) Suppose that S cq (<px) is a singleton and denote vx — v Vx . Then for all / € C(M) and v £ S Vx , 

lim -S n f(x) = vx(f), 

n—^oo Ji 

for v-a.e. x G M. 

(4) Suppose that S cq (ipx) is a singleton for X small enough and that g(X 1 Y) is C 1 ' 2 in a neighborhood of '(0, 0). 
Then the transport coefficients are defined and satisfy the Onsager reciprocity relations. 

(5) Suppose that the assumptions of {A) hold and write v = vq and = $q^. Then v(<f>(fi) = 0. Suppose in 
addition that v(Jf>( k )<f>n^) = 0(n _1 )/or n — > oo and that for some e > 0, 

sup -| \ogv(e~ Sn( - Y '^)\ < oo. 

Y£D t ,n>0 n 

Then the Fluctuation-Dissipation Theorem holds. 



59 



(6) Suppose that Y H > 1") ;s differentiable for all Y and let v € iSp x . 77ien f/ze Large Deviation Principle 
holds for &x w.r.t. (M,<fix, v) with the concave rate function Ix{s) — su PyeR«(i^ • s + q(X,Y)) which 
satisfies the GGC- symmetry Ix{s) = X ■ s + Ix(—s). 

The proof of Proposition 1 10. 131 (1) is the same as the proof of Proposition ! 10. 121 (1). The proof of the remaining 
statements is the same as the proof of the corresponding statements in Proposition l5.15l 

Assuming "chaoticity" one can say more. For example if cf>x is an expansive homeomorphism with specification 
and tpx and <&x are regular, then the entropy maps are upper-semicontinuous, S C q((fx — Y ■ $x) is a singleton 
for all Y (Theorem \TUM and the map Y H> q(X, Y) is everywhere differentiable (Theorem ll0.2l) . If {M,4>x) 
is a topological + transitive Smale space and (fx, $x are Holder continuous, then the map Y n- q(X, Y) is real 
analytic. In addition, Theorem l 1 0. 1 01 yields : 

Proposition 10.14 Suppose that (M,cf>x) is a topologically mixing Stnale space and that ipx, <&x> are Holder 
continuous for X in a neighborhood oj "0. Suppose also that g(X, Y) is C 1,2 in a neighborhood of (0, 0). Then the 
Fluctuation-Dissipation Theorem holds. 

10.4 Reference measure and physical equivalence 

In this section we investigate the relation between the observable a v introduced in the previous section and the en- 
tropy production observable a defined by equ. ( 13.151 ). Throughout the section we make the following assumptions: 

The homeomorphism <fi admit a time reversal § and ip £ Cr(M) is a potential. 

The metric d#(x,y) = d(i3(x),'0(y)) is equivalent to d, i.e., that there exists a constant C > such that 
C _1 d(x,y) < d#(x,y) < C d(x, y) for all x, y £ M. This requirement is a mild regularity assumption on 
It follows that d is also equivalent to d + dg, so that we can as well assume that d is isometric. 

u € S v is a ITRII reference state and the dynamical system (M,<fi,Lu) satisfies Assumption |(C)| with an entropy 
production observable a — £ Ul \ u o <fi <= Cr(M). 

The following key proposition relates a to u v . 



Proposition 10.15 Under the above assumptions one has a ~ a 



<p- 



Proof. Using the elementary identity B n (-d o <fi n 1 (x), e) = 'do cj) n 1 (B.„ (x. e)) and the fact that u> is lTRll we can 
write 

u(B n (d o <p n -\x),e)) = u- n+1 (B n (x,e)) = lo (c- s "- ia l Bn(x ^ . 
We derive the inequalities 

e s n -Mx)-(n-i)s c ^ Bn ^ e ^ < UJ {B n {do<p n - 1 {x) 1 e)) < e- s "- ia(x)+( "- 1)5 ^(B„(x,e)), 

where S e — sup x yeM M Xty s <e \cr(x) — a(y)\. With c = max Ig M l cr ( a; )l' we thus obtain, for arbitrary n e N, 
x £ M and e > 0, 

e _ nSr _ c < u (B n (0o<F-Hx),e)) eSna{x} < enSr+c _ (1Q g3) 

io{B n (x, e)) 

Since ui E S v and (S n <p)($ o (j) n ~ x {x)) = S n (<p o #)(&), the estimate (I10.73l l leads, for small enough e > 0, to 

Cn(e)- 1 < u:{B n {do4 > n ~ 1 {x),e))e- s ^ a ^ < C„(e). 
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Using again ( 110.731 ), we obtain 



which, combined with 010.831 , yields 

1 c 2 

-\S n {a-a v ){x)\ <<5 e + - + -logC n (e). 
n n n 

It follows that 

limsup - sup \S n (a - a v )(x)\ < S e , 
and the proof is completed by noticing that 5 e I as e I 0. □ 

Corollary 10.16 (1) The ES-functional of the dynamical system (M,(f>,u) 

l3(i4 e(a) = lim - logo^e-" 5 ™ 17 ), 
exists ant/ satisfies the ES-symmetry e(l — a) = e(a). 

/« //ze remaining statements, we assume that <f> is expansive with specification and that the potential tp is 
regular. 

(2) The ES-functional e(a) is everywhere differentiable. 

(3) The ES Fluctuation Theorem hold. 

(4) The system (M, 0, w) /zas a unique NESS uj + and for any f € C(M), 

lim -S n f(x) = u+(f), 

n— foo 77, 

holds for oj-a.e. x € M. 

(5) 77ze GC Fluctuation Theorem hold. 

(6) 77;e principle of regular entropic fluctuations hold. 

(7) Tjf w+ /n non-atomic, the system is entropy producing ijfui+ =/= uj+ o <&. 

Proof. (1) Proposition 110.31 yields the existence and the relation e(a) = P(<p — aa). Proposition 110. 1 51 implies 
(p — aa ~ <p — aa v from which Proposition 1 1 0.71 (2) allow us to conclude that e(a) = e v (a). Thus, the ES- 
symmetry follows either from Proposition ! 10. 121 (2) or directly from Proposition ^ .41 (2). 

Invoking Theorem ll0.81 (2)-(6) are direct consequences of Proposition ! 10. 121 In particular lj + = v v € S v so that 
(7) is a direct application of Proposition ! 10.1 11 (6). □ 

Remark. Suppose that (M, <ft) is a topologically + transitive Smale space and that a v and a are Holder continuous. 
Applying Propositions 110.101 we conclude that the ES-functional e(a) is real analytic. Moreover, <r„ and a are 
homologous, i.e., 

a v = a + ho <f) — h, 

for some h € Cr(M) which is unique up to an additive constant and also Holder continuous if (M, 0) is regular. 

Invoking Proposition ! 1 0. 1 3l similar results can be obtained for systems (M, <px , ) , depending on control param- 
eters, with ITRll reference states ux <E S vx . Note that if $>x is a continuous flux relation for the corresponding 
entropy production observable ax, then ( 110. 82b holds. We leave the details to the reader. 
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10.5 Markov chains 

We shall illustrate the results of this section on the simple example of a Markov chain with finitely many states. 
The set of states is Cl = {1, • • • , I}, 

M = Cl z = {x = {xj) je z | Xj e Cl}, 

with the usual product topology and <p ■ M — > M is the left shift, <j)(x)j = Xj+i ■ M is metrizable and a convenient 
metric for our purposes is 

d{x,y) = \ k ^\ 

where A e]0, 1[ is fixed and k(x, y) = inf{|j| | Xj ^ yj}. If d(x, y) < 1, then x = y . It follows that <p is 
expansive and that any r e]0, 1 [ is an expansive constant. Setting 

[x,y] = (. . . ,y- 2 ,y-i,x ,xi,. . .). 

one easily shows that (M, 4>) is a regular, topologically mixing Smale space. 

Any function / : M — >• R which depends only on finitely many xj's is Holder continuous. The map $(x)j = x_j 
is an isometric time-reversal. 

Let P = \pij]i,jen with pij > 0, J2j Pij — 1. be a transition matrix. By the Perron-Frobenius theorem there is a 
unique probability vector p = [pjjign with p i > and J^i Pi — 1 sucn tnat P P = p. 

The (two-sided) Markov chain with transition matrix P is the invariant Borel probability measure i/e5; such that, 
for any cylinder, 

C = {x e M | x k = ji,x k+1 = j 2 , ■ . . ,x fe+n _i = j n }, 

we have 

The assumption that p^ > implies that v is mixing w.r.t. (f>. 
For the potential 

tp(x) = logp X0X1 , (10.84) 

and the corresponding observable 

a v (x) = (p(x) - (p o d(x) = log PxoXl . 

PxqX-% 

one computes 

= PiPij lo SPv> v (°v) = PiPijlog^. 

Since B n (x, X k ) = {y e M \ y - L = x t for — k < i < k + n - 1} for k > 0, it easily follows that for any < e < 1 
there is a constant C e such that, for all n > 0, 

C- 1 <v(B n (x,e))c- s ^ <C e . 

We conclude that 

The partition ^ = ({a; e M | xo = i})ien has diam(^) = A < 1. Hence, the Kolmogorov-Sinai entropy can be 
computed from h v (<p) = h v {<\>, £) and a simple calculation leads to 

h v {4>) = - PtPv lo &Pii- 

i,jeS2 
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The pressure e v (a) = P(ip — aa v ), a real analytic function of a, is most easily computed from Equ. ( 110. 77K 
where we have 

Z n & - aa v ) = Z n ((l - a)<p + a<p o 0) = £ (p££j£ a *i) • ■ ■ fc>^J = trPS, 

with the matrix P Q = [py (a)], Pij (a) = p i J a p'j i . Since Pij(a) > 0, the Perron-Frobenius theorem applies to P Q 
and consequently e v (a) is equal to the logarithm of its dominant eigenvalue. In particular one checks 

P(tp)=0 = u(ip)+h v (<f>), 

so that v is the unique equilibrium state for the potential ip. Note that P* = Pi_ a . The resulting identity tr P™ = 
trP"_ Q provides an alternative proof of the symmetry e v (a) — e v {\ — a). 

Proposition 110. 1 1 1 implies that v{a v ) > and v(a v ) = iff v o ■& = v. The latter condition is easily seen to 
equivalent to PiPij = PjPji for all i,j £ f2. In other words, v(a v ) = iff the Markov chain satisfies detailed 
balance. Note that in this case a v = g o <f> — g ^ 0, with g(x) = \ogp XQ p X0X _ l . 

Let q = [<7i]i<=fi be a probability vector and u> the state uniquely determined by 



u) {{x £ M | x k = j k ,k = -to, ...,n}) = q jo (p joh ■ ■■Pj„_ 1 j n ) {Pj Q3 -i ■ "Pj. 



A simple calculation shows that oj is lTRII Moreover, lo £ S v provided qi > for all i £ fl. Thus, Corollarv ll0.16l 
applies to the ITRII system [M, <f>, ui). Note in particular that its entropy production observable 

a(x) = log , 

QxxPxxXq 

is homologous to a v . Explicitly, a v — a = h o tp — h with h(x) = log(q xo p XQX1 ). The unique NESS of the 
system is oj + = v and the system is entropy producing iff v does not satisfy detailed balance. 



Suppose that the transition matrix Px — [pij(X)] depends on the control parameters X £ M. N , We assume that 
the functions X M- pij(X) are C 2 and that pij(X) > for all X. Denote by p(X) = [J^pf)] the corresponding 
equilibrium vector. Let ipx be the corresponding potential, vx its equilibrium state and set ax — fx — fx ° 
We assume that detailed balance holds for X = so that ctq ~ 0. 

For each X £ l w let q(X) = [qi{X)\ be a probability vector such that qi{X) > and assume that q(0) = p(0). 
Construct the ITRII state u>x as above and denote by ax the corresponding entropy production observable. The 
detailed balance condition at X = implies that ljq = vq and hence <7o = 0. 

Setting 

Jo Pi](uX) 

and $x(%) = Fx(xo, x%) — Fx(%i,%o) we obtain a flux relation, 

a x ~ o x ~ X ■ & x , 

such that the map X i-> <f>x G C(M) is differentiable and <I>x 3»x- Arguing as before, the assumption 

Pij(X) > and the relation 

Z n (<p x -Y-<f> x )=trF(X,Y) n , 

where F(X, Y) = \p ij (X)e- Y < Fx ^- Fx( ' j ' i ^} imply that q(X, Y) = P x (fx-Y ■ * x ) is the logarithm of 
the dominant eigenvalue of P(X, Y). The perturbation theory of isolated simple eigenvalue further implies that q 
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is C 12 in a neighborhood of (0, 0) and all the conclusions of Propositions [TO. 1 31 and [1 0. 141 hold. Finally, for the 
family (M, 4>x,ujx, &x), one shows that 

g(X,Y) = g + (X,Y) = g(X,Y). 

Remark. The above approach can be used to discuss entropic fluctuations of a wide range of stochastic processes 
with suitable modifications to accommodate general (non-compact) state space, continuous-time, as well as Gibbs 
measures rather than Markov measures. The entropic fluctuations of Markov chains were first discussed by Kur- 
chan (Kul '|, Lebowitz and Spohn |LS2|, and Maes et al. llMalllMRVIIMNIl who used the path measure approach 
and the Gibbsian formalism. 



11 Anosov diffeomorphisms 

Let M be a compact connected smooth Riemannian manifold with a given Riemannian metric and let u> be the 
induced volume measure on M. We denote by Diff fc (Af) the set of all C k diffeomorphisms of M equipped with 
the usual C fe -topology. 

4> e DiS 1 (M) is called Anosov if M is a hyperbolic set, i.e., if there exist constants < A < 1, K > 0, and a 
decomposition of the tangent bundle 

TM = E u ®E S , (11.85) 
into £>(/>-invariant unstable and stable subbundles, such that for each x € M and every n € N 

II D x cj) n \ Ea || < KX n , \\D x <}>- n \ FU \\<K\ n , (11.86) 

x x 

The above bounds should hold for some norm equivalent to the Riemannian metric of M and there always exists 
such a norm for which K = 1. 

The set Ak(M) of all C k -Anosov diffeomorphisms of M is an open subset of Diff fc (M) (which can be empty). 
An Anosov diffeomorphism is called transitive if for any two non-empty open sets U and V and any N > there 
exists n > N such that 4>~ n (U) n V ^ 0. Any Anosov diffeomorphism on a torus T™ is transitive. More generally, 
it is conjectured that all Anosov diffeomorphisms are transitive (see [KH] for various partial results). 

In this section we consider dynamical systems (M, 4>, ui) where <fi g A2(M). We use freely notations and results 
from Section [Tol 

Let D x (f> : T X M — > T^ X )M be the derivative map of (f> at x. The entropy production observable of (M, <f>, ui) is 

<j(x) = -\ogD(x), 

where D(x) — \detD x <t>\ is the Jacobian of <fi at x. Setting D u (x) = \ detD x (j>\ Eu |, D s (x) — \ detD x (t>\ EB |, we 
shall consider the potential 

ip(x) = -logD u (x), 
which is known to be Holder continuous for <fi € A2(M) (see [Bo2 PS1 1). 

We shall say that the system (M, (j),uj) is ITRIl if it satisfies the conditions of Section [23] with a continuous time 
reversal •&. 

Remark 1. If i? G Diff 1 (M) then D x $ provides an isomorphism between E x u and E^^E^^y in particular the 
stable and unstable subbundles have the same dimension. Moreover, 

log D u otf = -log D s of 1 , (11.87) 
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so that a is homologous to a v = (p — (p o 

a 9 -a = logD s - logZ) s o (jr 1 . 
We stress however that we shall not assume § to be of class C 1 in the following. 

Remark 2. One can always construct lTRII Anosov systems starting with an Anosov system (M, <f>, uj) and applying 
the construction described at the end of Section l2~5l The time reversal obtained in this way is C°°. 

The following classical result is known as the Volume Lemma ([Bo2 KH|): 

Theorem 11.1 If '<f> £ A^{M) then u> £ S v . More precisely, for sufficiently small e > there exists C' e > such 
that, for all x £ M and n > 0, 

C- 1 <u>(B n (x,e))e- s ^W < Ce . 

Note that, by Corollary [T061 this implies P(tp) = 0. Set 

e„(a)= W (e^ 5 " CT ). 
If (M, (j>, u>) is lTRII the finite time ES-theorem (Proposition l3.5l ) yields that 

e„(a) = e„(l - a). (11.88) 
Proposition 110.31 immediately imply the existence of the ES-functional. 

Proposition 11.2 Suppose that 4> G Ai(M). Then for all a £ M, 

e(a) = lim — log e n (a) — P{<p — aa). 

n— >ao fi 

If(M, <p,u) isUMthen 

e(a) = e(l-a). (11.89) 



Remark. The symmetry dl 1.891 ) is forced by the finite time symmetry ( II 1.881 l and so Relation ( II 1.87b (and the 
fact that -d is C 1 ) is not used. One gets a direct proof of dl 1.891 ) based on (II 1.87b using the variational principle for 
the pressure in the same way as in the proof of Proposition ! 10. 121 (2). 

To improve Proposition 11 1.21 we need to assume more. We recall another classical result in the theory of Anosov 
diffeomorphisms: 

Theorem 11.3 Suppose that <f> € A%(M) is transitive. Then (M, <p) is a topologically mixing Smale space. 

In particular <f> is expansive with specification and Theorem ll0.8l together with the fact that P(<p) — 0, yield that 
4> has a unique equilibrium state v v G S v . Furthermore, by Corollary 1 10.41 

lim -S n f(x) = v v {f), 

n— foo Ji 

for all / £ C(M) and w-a.e. x. Thus, u> + = v v is the unique NESS of (M,<fi,uj). Proposition 110.31 yields that 
GC-functional exists and is equal to the ES-functional. 
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Proposition 11.4 Suppose that <fi S A2 (M) is transitive. Then for all a € M, 



e+(a) = lim — logw + (e 



-aS„a 



) = P{if-aa). 



n— foo 72 



/n particular, e(a) = e+(a) and (M, 0, /zas regular entropic fluctuations. 

The analyticity of e(a) and e+ (a) (Theorem llO.101 ) yields the respective Large Deviation Principles for the entropy 
production observable. 

Regarding the strict positivity of entropy production we have the following result. 
Proposition 11.5 (1) Suppose that 4> G A2{M) is transitive. Then oj+{<j) = ijfoJ+ <C u). 

(2) Suppose that all Anosov diffeomorphisms of M are transitive (for example, M — T n ). Then there is an open 
dense set A C Ai{M) such that for all (j) G A, w + ((t) > 0. 

Proof. (1) If w+ <C ui, then Corollarv l5.5l implies that oj+(u) = 0. The other direction follows from the result of 
Ruelle IIRu61 . (2) follows from (1) and the stability result of Sinai J2D which states that for an open dense set of 
(j)'s in A 2 (M) the NESS is singular w.r.t. uj. □ 

The resonance interpretation of e+(a) follows from Theorem 1 1 0. 1 01 (6). The resonances interpretation of e(a) 
follows from recent results of Baladi and Tsujii | Ba2 Ba3 ] and Gouezel and Liverani BBKLI IGLU IGL2I ILi2l ILTl 
on the spectrum of transfer operators in anisotropic Banach spaces and on the zeta function for Anosov maps. 

We now turn to the discussion of linear response theory for Anosov diffeomorphisms. Let k > 4 and let X <f>x 
be a C k map from some neighborhood of the origin in M. N into Ak(M) such that cj is an invariant state for 
4> = (f>o (Ak{M) is a Banach manifold so the notion differentiability makes sense). The map X H> ax = 
— log \detDtpx I S Cr(M) is C k ~ x . We shall consider only the flux relation 



Theorem 11.6 Suppose that all Anosov diffeomorphisms of M are transitive. Let k > 4 and let X t— > 4>x G 
Ak(M) be a C k map from some neighborhood of the origin in such that uj is an invariant state for cf> ~ 0o- 
Suppose that (M, 4>x,itj) is \TRl\ with a time-reversal independent of X . Then the Fluctuation-Dissipation Theorem 
holds: the transport coefficients 




Clearly, X ^ e C R (M) is G 



'fc-2 



are defined and satisfy the Onsager reciprocity relations 



Ljk — Lkj ■ 



For some A, B > and all n, \u($W$%>)\ < Ae~ B ^ and the Green-Kubo formula 



holds. The Central Limit Theorem holds for $ with covariance matrix [Djk] = 2[Ljk]. 



66 



Proof. (A/, <j>x) is a topologically mixing Smale space for X small enough. Arguing as in Propositions 111.21 and 
II 1.4| we deduce that 

g(X, Y) = g+ (X, Y) = P x (<p x -Y ■ * x ) . 

Ruelle [Ru3| has proven that the map X H> <px € C a (M) is C fe ~ 2 . Combining this result with Theorem 7 in 
HKKPW1 one deduces that (X, Y) H> g(X, Y) is C 1,2 in the neighborhood of the origin and the result follows 
from Proposition ll0.14l □ 

We finish with some remarks. 

Remark 1. The Green-Kubo formula and Onsager reciprocity relations for Anosov diffeomorphisms were first 
proven in [GR|. This proof was based on explicit computations and the differentiation formula established in 
llRulH . 

Remark 2. The proof of Theorem 11 1.61 looks deceptively simple. It stands on the shoulders of deep results 
established in llRu3l IKKPWI . 

Remark 3. The linear response theory for Anosov diffeomorphisms can be also established starting with the finite 
time Green-Kubo formula and following the strategy outlined in Section [5721 Obviously, the assumptions of the 
finite time linear response theory discussed in Section l4~3l hold under the conditions of Theorem ll 1.61 Under the 
same conditions Ruelle !IRu3H has proven that the functions X n- lox+ P ) are differentiable (see also MKKPWI 
and 1GL2ID . That the limit and derivative in the expression ( 15.41b can be interchanged follows from the results of 
Gouezel and Liverani MGL2I lLi3l (see also |Po|). 

Remark 4. Most of the results in this section extend to a certain class of Anosov flows for which sufficiently fast 
mixing has been proved, such as contact Anosov flows and flows with smooth stable and unstable foliations, see 
HBGMI IGel iDoTl lDo2l ILi 1IIBLL 



67 



Table of Abbreviations and Symbols 



(Q 


Equivalence of w and w< 


pas 


se|31 


(El) 


Regularity of 1 1-> A, , j, , 


pas 


je|I3 


(E2) 


Boundedness of a 


page [14] 


(E3) 


Existence of (er)+ 


pas 


ieM 


(M) 


Group property of (fr 


P a t 




(F2) 


Measurability of (t, x) i— > (x) 


pas 




(NESS1) 


M is a complete separable metric space 


pat 


?e|37] 


(NESS2) 


Continuity of (t, x) h-> 0*(x) 


pas 




(NESS3) 


Continuity of a 


pas 
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pas 
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pas 
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(T4) 


Existence of + 
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(T5) 


= Oft" 1 ) 

\ t / v j 
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ES -functional 


Evans-Searles functional e(a) 


nac 

r u ; 


;e[3Ql 


GC-functional 


Gallavoti-Cohen functional e+ (a) 


pas 


;e|38] 


GES -functional 


Generalized Evans-Searles functional g(X, Y) 


pas 


reEU 


GGC-functional 


Generalized Gallavoti-Cohen functional g+ (X, Y) 


pas 


je@3 


NESS 


Non-Equilibrium Steady State(s) 


pas 


?e|37] 


TRI 


Time-Reversal Invariance 


pas 


;ei 
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r 1 f*w(/.)da 


pas 
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(/> + 
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pas 
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Entropy cocycle 


pas 
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pas 
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Li = [Ljkt] 


Finite time Onsager matrix, transport coefficients 


pas 
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Onsager matrix, transport coefficients 


pas 
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Normal states w.r.t. uj 


pas 
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Invariant states 
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Radon-Nikodym derivative 
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Mean entropy production rate 
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